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FOREWORD 

This  report  describes  an  analytical  and  experimental  study  of  three  pavement 
configurations:  anchored,  elevated  anchored,  and  conventional  (rectangular 
slab).  The  purpose  of  the  analytical  study  was  to  present  thermal,  mechanical, 
and  thermo-mechanical  properties  of  materials  in  a  form  suitable  for  computer 
programming,  and  to  describe  environmental  and  mechanical  properties  of  the 
anchored  and  the  conventional  pavement  configurations.  Three  types  of 
analysis  were  conducted:  heat  transfer,  thermal  stress,  and  mechanical 
stress.  The  mechanical  stress  was  evaluated  by  two  methods:  finite  elements 
and  beam  theory.  These  analytical  results  are  compared  with  experimental 
results.  The  heat  transfer  analysis  was  used  to  evaluate  temperature  distri- 
butions within  the  pavement  and  the  soil  to  refine  the  mechanical  stress 
analysis. 

This  report  is  one  of  several  resulting  from  a  research  contract,  "New 
Structural  Systems  for  Zero-Maintenance  Pavements,"  issued  to  Dames  &  Moore 
by  the  Office  of  Research  and  Development  of  the  Federal  Highway  Administration, 
The  objective  of  this  research  study  is  to  identify  and  assess  the  potential 
of  new  and  innovative  structural  concepts  and  systems  to  serve  as  "Zero- 
Maintenance"  pavements.  An  interim  report,  "Unique  Concepts  and  Systems  for 
Zero  Maintenance  Pavements,"  Report  No.  FHWA-RD-77-76,  provides  an  updated 
state-of-the-art  and  comprehensive  review  of  each  of  the  three  major  struc- 
tural components  of  a  pavement  system:  the  subgrade,  the  base  and  subbase, 
and  the  pavement  surface.  Three  additional  manuals  are  being  prepared  to 
complete  the  study:  a  user  manual  for  the  computer  program,  a  construction 
specifications  manual,  and  an  analysis  manual.  These  manuals  are  tentatively 
scheduled  for  release  and  publication  by  the  FHWA  in  early  1980. 

Copies  of  the  report  are  being  distributed  by  transmittal  memorandum 
primarily  to  research  and  development  audiences. 

Charles  F.  Sch^P^ey 
Director,  Office  of  Research 
Federal  Highway  Administration 
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i6.  Abstract    New  Structural  Systems  for  Zero-Maintenance  Pavements^-The.  purpoge  of  this 
study  is  to  investigate  the  feasibility  of  designing  and  jconstructfrig1  cost-tjffective 
"Zero-Maintenance"  highways 

Volume  1:  Analytical  and  Experimental  Studies  of  an  Anchored  Pavement.  This 
report  documents  an  investigation  of  the  design  feasibility  and  construction  cost- 
effectiveness  of  an  anchored  pavement  concept  for  Zero-Maintenance  highways.  An 
analytical  model  is  designed  to  verify  computer  program  results  and  to  investigate 
construction  methods  for  a  full-scale  highway  section.  The  purpose  of  the  analytical 
study  is  (1)  to  present  thermal,  mechanical,  and  thermomechanical  properties  of 
typical  materials  in  a  form  easily  adaptable  to  computer  programs,  and  (2)  to  describe 
environmental  and  mechanical  properties  of  a  conventional  slab  and  an  anchored  pave- 
ment in  both  continuous  and  jointed  configurations.  The  two  pavements  were  subjected 
to  heat  transfer,  thermal  stress,  and  mechanical  stress  analyses.  The  anchored  pave- 
ment offers  two  distinct  advantages  over  a  conventional  pavement  -  deflections  are 
lower  and  more  uniform,  and  stresses  in  the  soil  are  lower  and  distributed  more  widely 
by  the  rigid  anchors.  Subgrade-related  failure  is  less  likely  to  occur  if  loads  are 
transmitted  deeper  within  the  subgrade.  Three-dimensional  finite  element  analysis  is 
considered  to  be  the  most  efficient  technique  for  examining  the  significance  of 
environmentally  induced  stress.  The  use  of  the  finite  element  method  is  anticipated 
as  more  advanced  analytic  techniques  are  developed. 

This  volume  is  the  first  in  a  series.  The  others  in  the  series  are:  Volume  2: 
Analysis  of  Anchored  Pavements  Using  ANSYS  and  Volume  3:  Anchored  Pavement  System 
Designed  for  Edens  Expressway.  Abstracts  of  these  volumes  are  included  on  page  ii  of 
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Volume  2:  Analysis  of  Anchored  Pavements  Using  ANSYS 

This  report  is  a  manual  which  provides  a  set  of  procedures  to  evaluate 
the  response  of  an  anchored  pavement  subjected  to  vehicle  static  loads, 
moisture  variation  in  the  subgrade,  and/or  temperature  variation  through 
the  surface  of  the  pavement.  These  procedures  include  two  computer  pro- 
grams known  as  FEMESH  and  ANSYS.  The  FEMESH  program  generates  rectangular 
meshes  in  either  a  two  or  three  dimensional  coordinate  system  for  any 
prespecified  number  and  spacing  of  nodes.  The  ANSYS  program  evaluates  the 
stresses,  strains,  and  deflections  at  all  elements  in  each  material  in- 
cluded in  the  analytical  model.  The  program  can  be  used  for  any  number  of 
different  materials  in  any  direction.  In  the  analysis  of  heat  transfer, 
the  program  provides  the  distribution  of  temperature  as  a  function  of  time 
at  predesignated  points.  The  program  is  versatile  and  capable  of  solving 
complex  geometrical  structures  supported  on  a  geologically  complex  earth 
mass.  The  behavior  of  an  anchored  pavement  section  is  evaluated  with  sets 
of  computer  programmed  mechanistic  models.  The  manual  was  written  to 
minimize  reference  to  other  publications. 


Volume  3:  Anchored  Pavement  System  Designed  for  Edens  Expressway 

This  report  provides  an  analysis  example  of  an  actual  pavement  and  the 
cost  estimate  using  the  anchored  system.  The  actual  pavement  is  the  Edens 
Expressway  in  Chicago.  The  report  provides  the  response  of  the  Edens 
Expressway  subjected  to  mechanical  and  environmental  loads  using  the  an- 
chored pavement  concept.  The  mechanical  and  thermal  properties  of  materi- 
als that  could  be  encountered  in  future  reconstruction  of  Edens  Expressway 
are  presented  in  a  consistent  form  for  computer  programming.  These  proper- 
ties are  viewed  as  typical  design  values  during  investigation  of  pavement 
response.  The  behavior  of  the  anchored  pavement  under  induced  temperature 
loads  and  weakening  of  subgrade  (by  thawing  action)  is  clearly  demon- 
strated. This  report  will  enable  application  of  the  anchored  pavement 
concept  by  any  road  with  heavy  traffic.  The  example  problem  provides  the 
input  parameters  of  materials  and  loads  for  the  analysis,  the  generation  of 
finite  element  mesh,  and  the  results  of  the  analysis.  The  computer  program 
ANSYS  was  used  for  this  study  (the  manual  for  the  use  of  the  program  is 
presented  in  Vol.  2  of  this  series  of  reports). 
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CHAPTER  I- 
INTRODUCTION 

Serious  concern  over  the  costs  of  highway  systems  has  raised  a 
question  as  to  the  feasibility  of  designing  and  constructing  cost  effective 
"zero  maintenance"  highways.  For  this  project,  the  Federal  Highway  Adminis- 
tration (FHWA)  has  defined  zero  maintenance  as  that  requiring  no  structural 
work  for  20  years,  then  only  minor  work  for  the  next  10  years. 

As  a  result  of  research  sponsored  by  FHWA,  several  structural  con- 
cepts have  been  proposed  (Saxena,  Wang,  Udwari,  and  Rosenkranz,  1977)  for 
zero  maintenance  performance.  These  include  pile  support,  edge  stiffening, 
thick  cellular  systems,  waffle-type  systems,  modified  conventional  systems, 
and  a  flexible,  "floating,"  V-shaped  pavement. 

A  limited  study  of  various  concepts  showed  an  anchored  pavement 
(Figure  1)  to  be  promising  because  it  uses  a  similar  amount  of  structural 
material  as  present  systems,  does  not  pose  great  construction  difficulties, 
and  probably  requires  little  subgrade  preparation. 

1.1  SCOPE  OF  STUDY 

Investigation  of  the  proposed  anchored  system  was  divided  into 
two  sections:  Model  Study  and  Analytical  Study. 

The  scope  of  the  model  study  was  twofold:  first  and  foremost, 
to  verify  the  analytical  model  (computer  program)  results  and  second,  to 
investigate  construction  methods  and  problems  for  a  full-scale  highway  section. 
The  model  study  should  also  help  determine  specifications  for  a  full-scale 
test. 

The  analytical  study  had  two  intents:  the  first  was  to  present 
thermal,  mechanical,  and  thermomechanical  properties  of  typical  materials 
encountered  in  highway  construction' and  pavement  design  in  a  consistent  form, 
easily  adaptable  to  computer  programs.  These  properties  were  compiled  after 
a  comprehensive  literature  search  to  ascertain  typical  values  and  expected 
ranges.  General  trends  as  well  as  actual  values  are  discussed  in  this  report. 
The  second  was  environmental  and  mechanical  analyses  of  a  conventional  slab 
and  an  anchored  pavement  (Figure  1)  in  both  continuous  and  jointed  configurations 
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Figure  1.     Configuration  of  conventional  and  anchored  pavements  studied. 


In  the  heat  transfer  analysis,  extreme  temperatures  were  simulated 
to  determine  temperature  distributions  within  the  structural  elements  and 
supporting  subgrade;  results  were  compared  with  field  test  data  reported 
in  the  literature.  Thermal  stress  analysis,  based  on  temperature  distributions 
derived  in  the  heat  transfer  analysis,  was  also  conducted  and  results  were 
compared  with  those  reported  in  the  literature.  To  complete  the  environmental 
analysis,  the  effects  of  moisture  and  its  transport  through  the  systems  were 
studied. 

Pavement  response  to  mechanical  loads  imposed  by  vehicles  was  investi- 
gated for  all  configurations.  The  results  from  the  mechanical  stress  analysis 
were  superimposed  on  the  results  from  the  environmental  stress  analysis  so 
that  the  relative  magnitudes  of  stresses  and  displacements  could  be  seen. 
The  most  detrimental  load  (worst  case)  in  each  system  is  presented  so  that 
the  gross  behavior  of  each  system  can  be  understood.  Finally,  the  relative 
advantages  and  disadvantages  of  each  system  are  discussed. 

Two-dimensional  finite  elements  were  used  for  heat  transfer  and 
stress  analysis,  and  three-dimensional  finite  elements  for  mechanical  stress 
analysis.  A  computer  program  entitled  Engineering  Analysis  System  (ANSYS) 
developed  by  Swanson  (1975)  was  used  because  of  its  versatility.  With  ANSYS, 
several  types  of  analysis  are  possible  including  static,  dynamic,  and  heat 
transfer  of  two-  and  three-dimensional  structural  systems.  A  large  element 
library  is  also  available,  and  material  and  geometric  nonlinearities  are 
possible. 

The  finite  element  method  was  chosen  because  it  lends  itself  to 
analysis  of  complex  nonlinear  and  transient  field  problems  in  several  dimen- 
sions. Such  analysis  would  be  otherwise  impossible  without  costly  experimental 
investigations. 

The  investigation  is  reported  in  the  following  sequence: 

Chapter  II  -  Investigations.  Presents  experimental  study 
of  the  anchored  pavement  62  in.  (157.5  cm) 
long.  For  comparison  purposes,  a  conventional 
slab  was  tested  under  similar  conditions. 
The  results  of  a  series  of  other  experiments 
are  also  presented. 


Chaptar  III  -  Verification  of  Computer  Program.  Presents 
the  types  of  verification  and  the  numerical 
values  derived  from  experiments  and  actual 
field  tests  to  form  a  comparison  for  the  analyti- 
cal solution. 

Chapter  IV  -  Analysis  of  Pavement  Concepts.  Presents  material 
characterizations,  analysis  of  environmental 
loads,  analysis  of  mechanical  loads,  and  superpo- 
sition of  mechanical  and  environmental  effects. 

Chapter  V   -  Conclusions  and  Recommendations. 

The  report  also  has  a  bibliography  and  three  appendices  dealing 
with  soil  properties  in  experimental  work  and  finite  element  formulations. 

1.2  RELATED  DOCUMENTS 

This  document  is  developed  so  that  it  can  be  used  with  a  minimal  reference 
to  other  materials.  Two  subsequent  reports  which  include  information  on  anchored 
pavement  system  concepts  are  Volume  2:  Analysis  of  Anchored  Pavements  Using 
ANSYS  and  Volume  3:  Anchored  Pavement  System  Designed  for  Edens  Expressway.  See 

page  ii  of  this  volume  for  information  on  the  limited  availability  of  these  sub- 
sequent volumes. 


CHAPTER  II 
EXPERIMENTAL  INVESTIGATIONS 

2.1  GENERAL 

The  response  of  subgrade  to  load  applied  on  conventional  pavement 
slabs  has  been  previously  investigated.  Experimental  investigations  of  the 
structural  behavior  of  rigid  pavements  were  made  at  Arlington  Experimental  Farm, 
Va.  (Teller  and  Sutherland,  1935,  1936,  1943)  and  at  Iowa  State  Engineering 
Experimental  Station  (Spangler,  1935,  1937).  These  experiments  and  the  work  by 
the  U.S.  Army  Corps  of  Engineers  resulted  in  modification  of  Westergaard's 
formulas  (1927)  for  the  most  critical  conditions  when  the  slab  corners  are  warped 
upwards. 

Full-scale  tests  were  performed  at  Schiphol  Airport  in  Holland  (Vander 
Veen,  1953)  and  on  Interstate  80,  near  Ottawa,  Illinois,  as  part  of  the  AASHO 
Road  Test  Program.  The  latter  generated  many  studies.  Because  the  above 
mentioned  were  tests  conducted  in  the  field,  it  was  difficult  to  separate 
moisture  and  temperature  effects  from  load  effects.  Model  tests  were  performed 
under  controlled  conditions  to  study  the  effect  of  stresses  due  to  load  alone  by 
Vesic  and  Saxena  (1970).  In  these  tests,  behaviors  of  slabs  resting  on  an 
elastic  subgrade  of  known  properties  were  relatively  close  to  those  predicted  in 
analytical  studies. 

2.2  MODEL  DESCRIPTION  AND  INSTRUMENTATION 

The  model  tests  for  this  investigation  were  l/20th  scale  and  involved 
anchored  and  conventional  aluminum  slabs  of  similar  lengths;  a  subgrade  of  known 
properties;  a  container  tank  for  the  soil  and  for  conducting  the  experiment;  and 
loading  and  measuring  equipment. 

2.2.1  Slabs 

Both  the  anchored  and  conventional  slabs  were  62  in.  (157.5  cm)  long, 
21.63  in.  (55  cm)  wide,  and  0.5  in.  (1.27  cm)  deep  (Figure  2).  For  the  anchored 
slab,  there  was  an  anchor  near  each  edge  and  a  central  anchor.  The  two  edge 
anchors  measured  1.88  in.  (4.8  cm)  deep  and  0.625  in.  (1.6  cm)  wide.  The  central 
anchor  was  also  1.88  in.  (4.8  cm)  deep,  but  was  1.25  in.  (3.2  cm)  wide.  The 
lengths  of  the  anchors  and  the  slab  were  obviously  the  same.  The  anchors  were 
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attached  to  the  slab  by  screws  centered  at  5  in.  (12.7  ,n).     The  top  and  bottom 
views  of  the  slab  are  shown  in  Figure  3. 

This  series  of  tests  were  conducted  to  keep  the  slab  5  in.  from  the 
edges  of  the  container  in  which  the  experiment  was  being  performed  (Figure  3). 
Tests  on  conventional  slab  were  performed  using  the  62  in.  (157.5  cm)  slab 
without  the  anchors. 

2.2.2  Subgrade 

The  subgrade  used  was  42%  kaolinite  clay,  42%  silica  sand,   and  16% 
water  (by  weight).     The  subgrade  was  classified  as  a  silty  clay  with  a  plasticity 
index  of  16  and  an  optimum  water  content  of  about  8%.     The  soil  was  mixed  with  an 
over-the-optimum  water  content. 

2.2.3  Container 

A  steel  tank  was  used  for  the  tests.  The  tank  was  divided  into  two 
areas,  testing  and  sampling  (Figure  4).  The  testing  area  *as  73  in.  (185.4  cm) 
long,  36.25  in.  (92.1  cm)  wide,  and  22.75  in.  (57.8  cm)  deep.  The  sampling  area 
was  27.88  in.  (70.8  cm)  long,  with  other  dimensions  being  the  same. 

A  soil  mixer  and  a  compactor  were  used  to  prepare  the  subgrade.  The 
silty  clay  subgrade  was  mixed  in  100-1 b.  (45.36  kg)  batches  then  deposited  in  the 
tank.  The  subgrade  was  compacted  in  2-in.  (5.08  cm)  layers.  For  a  uniform 
compaction,  two  passes  of  the  compactor  per  layer  were  found  to  be  sufficient. 
Before  placing  another  layer,  the  top  0.5  in.  (1.27  cm)  of  the  previous  layer  was 
raked  to  ensure  proper  bond.  Every  effort  was  made  to  have  uniform  compaction  in 
all  layers.  The  leveling  of  the  upper  layer  was  done  with  precision  to  ensure 
proper  contact  between  the  slab  and  subgrade. 

Anchors  were  placed  before  the  leveling  of  the  upper  layer,  then  the 
slab  was  screwed  to  them.  To  avoid  disturbing  the  soil,  this  method  was  con- 
sidered better  than  welding  the  anchors  to  the  plate.  Crossbars  were  utilized  to 
keep  the  anchors  in  position  during  compaction.  The  crossbars  were  later 
removed. 

2.2.4  Loading  and  Measuring  Equipment 

The  loading  platform  (Figure  5)  was  designed  to  represent  a  l/20th 
scale  model  of  two  rear  axle  trucks  with  four  tires  per  axle  and  a  maximum 
capacity  of  18,000  lb.  (8165  kg)  per  axle.  The  18,000-lb.  axle  load  on  a  l/20th 


a)  Top  view. 


b)  Bottom  view. 


Figure  3.  Views  of  a  62-in.  (157.5-cm)  long  anchored  slab, 
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scale  model  was  equal  to  a  scale  45-lb.  (20.4  kg)  load;  however,  greater  loads 
were  used  within  the  elastic  range  to  investigate  the  response  of  the  slab- 
subgrade  system.  For  load  application,  a  manual  jack  with  an  attached  load  cell 
was  used.  A  strain  indicator  unit  and  balance  unit  were  used  to  measure  the 
strains  of  10  foil  strain  gages  fixed  at  various  locations  on  the  slab.  The  foil 
gages  were  used  in  a  quarter  bridge  circuit  with  an  internal  120  ohm  dummy. 

The  following  notations  were  used  for  the  series: 

1.  Anchored  slab,  62  in.  (157.5  cm)  long,  on 
soil  subgrade  —  AS62 

2.  Conventional  slab,  62  in.  (157.5  cm)  long 
on  soil  subgrade  —  S62 

3.  Anchored  slab,  62  in.  (157.5  cm)  long  elevated  by 
0.5  in.  (1.27  cm)  from  soil  subgrade  —  ASE62 

2.3  TEST  DETAILS 

2.3.1  Anchored  Slab  on  Soil  Subgrade  -  AS62 

The  soil  was  placed  using  the  technique  described  in  Section  2.2. 
After  the  final  layer  was  placed  and  leveled,  a  film  of  oil  was  spread  over  the 
area  not  covered  by  the  slab,  to  prevent  evaporation  of  moisture.  This  area  was 
then  covered  with  thin  plastic  sheets. 

A  cross  section  and  a  plan  view  of  the  slab  on  the  subgrade  soil  are 
shown  in  Figures  6  and  7  respectively.  The  location  of  the  dial  gages  for 
measuring  surface  deflections  of  the  anchored  slab  are  shown  in  Figure  8.  The  16 
numbered  points  where  the  deflection  was  measured  were  used  to  draw  transverse 
and  longitudinal  deflections  along  sections  Tl-Rl,  T2-R2,  and  T3-R3,  and  along 
sections  Ll-Nl,  L2-N2,  and  L3-N3  respectively,  also  shown  in  Figure  8. 

Ten  foil  gages  were  placed  on  the  slab  to  measure  strains  and  deduce 
the  stresses  and  bending  moments.  Five  foil  gages  were  placed  on  top  surface  of 
the  slab,  two  on  the  bottom  surface,  and  the  other  three  on  the  anchors 
(Figure  9). 

The  aluminum  material  in  the  experimental  range  were  linear,  iso- 
tropic, and  elastic.  Therefore,  two  independent  elastic  constants  (E  =  modulus 
of  elasticity  and  v  s  Poisson's  ratio)  were  necessary. 
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Figure  6.  Cross  section  of  AS62  on  subgrade 
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VARIOUS  POSITIONS  OF  LOAD  WITH  FIXED  DIAL  GAGES  ON 
S62,  AS62.  AND  ASE62 
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VARIOUS  POSITIONS  OF  LOAD  AND  LOCATION  OF  LINES 
WHERE  SURFACE  DEFLECTION  ARE  PLOTTED  FOR  S62,  AS62,  AND  ASE62 
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LINES  WHERE  SURFACE  DEFLECTION  HAVE  BEEN  PLOTTED 

V///X       LOAD  POSITION 
K7/4&r\V3  L0AD  P0S,T,0NS  1  AND  4  .  CROSS  HATCHED  AREAS 


o 

S62 
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ASE62 


BEING  COMMON 

DIAL  GAGES 

REPRESENTS  62"  LONG  SLAB  ON 
SOILSUBGRAOE 

REPRESENTS  62"  LONG  ANCHORED  SLAB  ON 
SOILSUBGRAOE 

REPRESENTS  62"  LONG  SLAB  ANCHORED  SLAB 
ELEVATED  0.5"  ON  SOIL  SUBGRADE 


Figure  8.    VARIOUS  POSITIONS  OF  LOAD,  DIAL  GAGES,  AND 

SURFACE  DEFLECTIONS 


14 


A-  7.38" -y 


NOTE: 

For  Locations  of  Oial  Gages  and  Lines  of 

Deflection,  See  Figure  3. 


Figure  9   .     POSITION  OF  FOIL  GAGES  ON  ANCHORED  SLAB 
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The  stresses  (Desai  and  Abel,  1972)  can  be  computed  as  follows: 


-  2- 


(1  +  v)(l  -  2v) 


1-v 


1-v 


l-\> 


The  bending  moments  (Popov,  1968)  can  be  also  computed  from  the  stresses  as: 


"Mx~ 

°z  lxy 

M 

y 

n     1 
c 

ax    yz 

M 

z 

fy  !zx 

where 


c  =  distance  of  foil  gages  from  the  neutral  axis 


The  modulus  of  elasticity  of  the  aluminum  used  was  10.5  x  10  psi 
2 
(738.150  kg/an  )  and  the  Poisson's  ratio  was  0.33. 

A  "line  preloading"  was  initially  used  at  the  center,  quarter  dis- 
tances, and  edges  of  the  anchored  slab.  The  preloading  was  considered  necessary 
to  prevent  loss  of  contact  between  the  anchored  slab  and  the  soil  and  to  bring 
the  soil  within  the  elastic  range.  The  preloading  is  shown  in  Figure  10.  The 
hysteresis  loops  (load  versus  deflection)  for  loads  at  the  center,  and  at  the  two 
edges  of  the  anchored  slab  are  shown  in  Figures  11,  12,  and  13  respectively. 

The  deflection  at  the  loading  point  was  measured  from  dial  gage  No.  1. 
After  preloading  and  ensuring  that  the  soil  was  in  the  elastic  range,  the  deflec- 
tion measuring  dial  gages  were  placed.  Dial  gages  Nos.  1  and  13  measured  0.0001 
in.  (0.003  mm)  per  division;  Nos.  11  and  15  measured  0.0005  in.  (0.013  mm)  per 
division,  and  the  others  measured  0.001  in.  (0.025  mm)  per  division. 
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Six  load  positions  were  usad  to  investigate  the  response  of  a  "contin- 
uous pavement"  and  of  pavement  at  a  joint  with  zero  transfer  load  (Figure  8). 
Figures  14,  15,  and  16  show  various  loading  configurations. 

Longitudinal  (Ll-Nl,  L2-N2,  L3-N3)  and  transverse  (Tl-Rl,  T2-R2, 
T3-R3)  surface  deflection  lines  for  load  at  position  1  and  position  3  were 
plotted  in  Figures  19  through  30. 

The  load  increment  was  250  lb.  (113.4  kg)  although  the  maximum  load 
applied  was  1500  lb.  (680.4  kg).  The  resulting  deflections  are  tabulated;  the 
plots  only  show  the  deflection  for  250-lb.  (113.4  kg),  500- lb.  (227  kg),  and 
750-1 b.  (340  kg)  loads.  The  measured  strains  are  shown  in  Appendix  A,  Tables 
5  through  16. 

2.3.2  Conventional  Slab  on  Subgrade  -  S62 

Soil  placement  and  compaction  procedures  were  similar  to  those  for 
AS62.  Great  effort  was  made  to  level  and  surface  the  top  layer.  The  slab  was 
then  placed,  and  the  remaining  soil  surface  was  covered  by  oil  and  plastic  sheets 
as  previously  described. 

Deflections  were  measured  at  16  locations,  as  shown  in  Figure  8.  Load 
positions  and  the  sections  where  deflections  are  plotted  are  shown  in  Figure  8. 
Six  foil  gages  were  attached  on  top  of  the  slab,  and  four  were  attached  on  the 
bottom  of  the  slab  (Figure  17). 

The  preloading  system  was  then  mounted  at  the  center  of  the  slab.  The 
dial  gage  No.  1  was  placed  to  measure  the  surface  deflection  beneath  the  load. 
The  maximum  load  applied  was  750  lb.  (340  kg)  with  increments  (or  decrements)  of 
250  lb.  (113.4  kg).  Figure  18  shows  the  load  versus  deflection  relation  for 
three  cycles  of  loadings  at  the  center  of  the  slab. 

Subsequent  loading  did  not  produce  any  appreciable  change  in  the 
hysteresis  loop.  Similar  hysteresis  loops  were  derived  for  line  preloading  at 
the  left  and  right  edges  of  the  slab. 

The  slab  subgrade  system  was  also  preloaded  at  the  left  and  right 
intermediate  sections.  After  the  preloading,  the  regular  loading  was  applied.  A 
maximum  load  of  750  lb.  (340  kg)  (load  increment,  250  lb.  =  113.4  kg)  was  used  at 
six  locations  (Figure  8).  Tables  II  through  16   (Appendix  A)  contain  measured 
strains  for  the  six  load  positions. 
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y.      A-.  7.38"-^ 


NOTE: 

Por  Locations  of  Oiai  Gages  anc 

Lines  of  Deflection,  See  Figure  8. 


Figure  17.      POSITION  OF  FOIL  GAGES  ON  CONVENTIONAL  SLAB 
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Figure  19.  Longitudinal  surface  deflection  along  line  Ll-Nl,  load  position  1, 

S62  -  AS62. 
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Figure  20.  Longitudinal  surface  deflection  along  line  L2-N2, 
load  position  1,  S62  -  AS62. 
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Figure  21.  Longitudinal  surface  deflection  along  line  L3-N3,  load  position  1 

S62  -  AS62. 
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Figure  22.  Transverse  surface  deflection  along  line  Tl-Rl, 
load  position  1 ,  S62  -  AS62. 
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Figure  23.  Transverse  surface  deflection  along  Line  T2-R2,  load  position  1, 

S62  -  AS62. 
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Figure  24.  Transverse  surface  deflection  along  line  T3-R3, 
load  position  1 ,  S62  -  AS62. 
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Figure  25.  Longitudinal  surface  deflection  along  line  Ll-Nl, 
load  position  3,  S62  -  AS62. 
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Figure  26.  Longitudinal  surface  deflection  along  line  L2-N2, 
load  position  3,  S62  -  AS62. 
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Figure  27.  Longitudinal  surface  deflection  along  line 
L3-N3,  load  position  3,  S62-AS62. 
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Figure  28.  Transverse  surface  deflection  along  line 
Tl-Rl,  load  position  3,  S62  -  AS62. 
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Figure  29.  Transverse  surface  deflection  along  line  T2-R2, 
load  position  3,  S62  -  AS62. 
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Figure  30.  Transverse  surface  deflection  along  line  T3-R3 
load  position  3,  S62  -  AS62. 
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^Figures  19,  20,  and  21  are  plots  of  the  longitudinal  surface  deflec- 
tions for  loads  at  position  1.  Figures  22,  23,  and  24  represent  plots  of  the 
transverse  surface  deflections  for  loads  at  positon  1.  Surface  deflection  for 
load  at  position  3  are  shown  in  Figures  25  through  30. 

The  above  plots,  being  most  representative  for  a  static  load  in  a 
conventional  pavement,  are  shown  in  this  report. 

2.3.3  Elevated  Anchored  Slab  -  ASE62 

The  procedure  followed  was  similar  to  that  for  the  AS62  and  S62  exper- 
iments. Before  the  slab  was  screwed  to  the  anchors,  an  0.5  in.  (1.27  cm)  space 
was  left  between  the  slab  and  the  subgrade.  Figure  31  shows  a  cross  section  of 
the  model,  and  Figure  32  shows  a  view  of  experimental  set  up. 

Sixteen  dial  gages  and  six  load  positions  were  used  to  study  the 
elastic  response  of  the  soil -structure  system  (Figure  8).  Location  of  sections 
where  surface  deflection  is  plotted  is  also  shown  in  Figure  8. 

Ten  foil  gages  as  used  in  Section  2.3.1  were  utilized  in  this  investi- 
gation (Figure  9).  Preloading  of  the  soil  was  performed  to  keep  the  soil  in  the 
elastic  range  during  experiment.  The  preloading  was  done  at  center,  at  quarter 
distance  on  the  left,  at  quarter  distance  on  the  right,  and  near  edge  on  left  and 
right.  The  preloadings  were  done  symmetrically  to  avoid  any  possible  lifting  of 
the  anchored  slab.  After  preloading,  the  loadings  platform  was  replaced  and  the 
slab  was  loaded  at  six  places,  as  described  previously.  The  maximum  load  applied 
was  300  lb.  (136  kg)  with  increments  of  100  lb.  (45.36  kg). 

For  position  1,  longitudinal  surface  deflections  are  plotted  in 
Figures  33,  34,  and  35,  and  transverse  surface  deflections  are  plotted  in 
Figures  36,  37,  and  38.  For  position  3,  longitudinal  and  transverse  deflection 
are  plotted  in  Figures  39  through  44. 

Strain  measurements  are  shown  in  Tables  17   through  22    Appendix  A. 

2.4  RESULTS 

As  shown  in  the  next  chapter,  the  trends  of  observed  deflections  of  the 
model  test  agree  well  with  the  analytical  results  for  the  anchored  slab.  In  the 
experimental  work,  the  conventional  slab  was  taken  as  a  basis  of  comparison. 

The  surface  deflections  of  the  anchored  slab  in  the  longitudinal 
direction  and  the  central  transverse  section  for  the  central  loading  position 

39 


(position  1)  are  about  one-third  of  those  obtained  for  a  conventional  slab.  At 
sections  T2-R2  (quarter  distance  from  edge)  and  T3-R3  (near  edge).  The  conven- 
tional slab  exhibited  significant  uplifting,  whereas  the  anchored  slab  had 
almost  no  uplift  or  insignificant  uplift  only  at  the  edge. 

Similar  trends,  that  is,  differential  magnitudes  of  the  order  of  one- 
third,  were  observed  for  the  edge  loading  (position  3)  as  shown  in  Figures  25, 
26,  27,  29,  and  30.  The  uplifting  of  the  center  was  very  pronounced  for  the 
conventional  slab  but  insignificant  for  the  anchored  slab  (Figure  28). 

The  anchored  slab  was  also  compared  with  the  elevated  slab  system.  The 
surface  deflections  beneath  the  load  for  position  1  of  the  anchored  slab  are 
about  one-third  to  two-thirds  of  the  deflections  for  the  elevated  slab 
(Figures  33  through  44).  No  significant  uplifting  was  observed  in  the  elevated 
slab. 

Figures  39  through  44  indicate  the  following: 

1.  For  a  load  of  250  lb.,  (113.4  kg)  the  deflections  of 

the  anchored  slab  in  full  contact  (AS62)  and  the  elevated 
anchored  slab  (ASE62)  were  found  to  be  very  similar. 
This  indicates  that  much  of  the  load  was  being  carried 
by  the  anchors  to  the  soil  beneath. 

2.  Higher  loads  on  the  elevated  slab  would  have  caused  punching 
failure  and  hence  could  not  be  used  in  the  experiment; 
however,  analytical  investigation  (though  not  performed 

in  this  report)  could  be  of  interest. 
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Figure  33.     Longitudinal  surface  deflection  along  line  Ll-Nl, 
load  position  1,  ASE62  -  AS62. 
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Figure  35.  Longitudinal  surface  deflection  along  line  L3-N3, 
"load  position  1,  ASE62  -  AS62. 
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Figure  36.     Transverse  surface  deflection  along  line  71-31, 
load  position   1,  ASE52  -  AS52. 
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Figure  37.     Transverse  surface  deflection  along  line  T2-R2 
load  position  1,  ASE62  -  AS62. 
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Figure  38.     Transverse  surface  deflection  along  line  T3-R3, 
load  position  1,  ASE62  -  AS62. 


48 


-0.005 


o 
i— * 

►— 


O 


0.005 


0.010 


0.015 


0.020 


0.025 


0.030 


_     1    in.   =  2.54  cm 
1   lb.   =  0.45  kg 


inAn 

ASF*? 

ask? 

(lbs.) 

100 

0 

200 

A 

?50 

• 

300 

a 

AS62  represents  62  in.   long 
anchored  slab  on  soil 
subgrade 

ASE52  represents  62  in.  long 
anchored  slab  elevated 


\ 


0.5  in.  on  soil  subgrad 


\ 


\ 


Figure  39.  Longitudinal  surface  deflection  along  line 
Ll-Nl,  load  position  3,  ASE62  -  AS62. 
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Figure  40.  Longitudinal  surface  deflection  along  line  L2-N2, 
load  position  3,  ASE52  -  AS62. 
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Figure  41.  Longitudinal  surface  deflection  along  line  L3-N3, 
load  position  3,  ASE62  -.  AS62. 
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Figure  42.  Transverse  surface  deflection  along  line  Tl-Rl, 
load  position  3,  ASE62  -  AS62, 
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Figure  43.     Transverse  surface  deflection  along  line  T2-R2, 
load  position  3,  ASE62  -  AS62. 
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Figure  44,  Transverse  surface  deflection  along  line  T3-R3, 
load  position  3,  ASE62  -  AS62. 
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CHAPTER  III 
VERIFICATION  OF  COMPUTER  PROGRAM 

3.1  INTRODUCTION 

Before  any  numerical  method  is  used  to  solve  complex  structural  or 
field  problems,  it  should  be  checked  against  similar  (though  often  simpler) 
situations  which  either  possess  closed-form  or  otherwise  known  solutions. 
Equally  acceptable  checks  can  be  made  against  experimentally  observed  results  or 
actual  field  conditions  provided  care  was  taken  to  ensure  accuracy  in  their 
measurement.  This  section  presents  the  types  of  verification,  mathematical 
modeling  techniques,  and  numerical  values  derived  from  actual  field  conditions 
and  experiments  used  to  form  a  comparison  for  the  solution  and  thereby  forecast 
the  relative  degree  of  accuracy  of  further  studies. 

The  verification  study  makes  use  of  the  heat  transfer  and  mechanical 

(static)  stress  analysis  options  available  in  the  computer  program  ANSYS.  As  a 

result,  test  programs  were  developed  specifically  to  lend  credence  to  these  types 
of  analysis. 

3.2  VERIFICATION  FOR  HEAT  TRANSFER 

A  literature  review  of  problems  associated  with  thermal  soil  mechanics 
produced  a  series  of  papers  describing  a  road  test  in  Norway  during  the  1957-69 
period  (Noss,  1973;  Nordal,  1973;  Frivik  et  al_. ,  1977).  Atmospheric  conditions 
were  monitored  and  temperature  profiles  within  the  ground  (tautochrones)  at 
various  hours  of  the  day  were  obtained.  Of  critical  importance  to  this  investi- 
gation was  the  additional  reporting  of  material  properties  -  thermodynamic  and 
mechanical  (Frivik  et  al_. ,  1977).  With  such  a  complete  set  of  data,  it  was 
possible  to  construct  a  mathematical  model  to  simulate  the  prototypical  road  test 
section  behavior. 

The  site  —  the  Vormsund  test  road  —  consists  of  several  sections  of 
pavements.  The  section  modeled  is  at  station  724  +  7.5;  the  profile  is  shown  in 
Figure  45.  The  temperature  profiles  were  taken  along  the  pavement  centerline,  to 
assure  a  near  one-dimensional  system.  Since,  however,  the  mathematical  model 
used  for  actual  analysis  would  be  two-dimensional,  it  was  decided  to  model  the 
test  road  by  two-dimensional  isoparametric  finite  elements,  the  same  types  used 
in  the  primary  study.  The  mesh  is  shown  in  Figure  45.  No  heat  flow  was  allowed 
on  vertical  boundaries;  the  lower  boundary  was  kept  at  a  constant  temperature. 
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These  boundary  conditions  simulate  a  soil  mass  of  infinite  extent.  The  lower 
boundary  condition  is  consistent  with  observed  values  (Figure  46)  and  the 
modeled  constant  temperature  at  the  100  in.  (254  cm)  depth  is  sufficient. 

The  input  surface  temperature  was  approximated  by  piecewise  linear 
segments  of  3  to  6  hours  (Figure  47).  The  analysis  is  nonlinear,  as  material 
properties  can  vary  with  temperature.  The  reported  values  of  thermal  con- 
ductivity and  specific  heat  vary  only  slightly  in  the  temperature  ranges  con- 
sidered, and  a  nominal  value  for  each  was  chosen.  The  input  values  are  given  in 
Table  1  and  are  considered  typical.  The  analysis  uses  a  time-stepping,  tentative 
procedure  that  updates  the  stiffness  (conductivity)  and  damping  (specific  heat) 
matrices  after  each  iteration.  The  analysis  assumed  a  uniform  temperature 
throughout  the  system  and  then  cycled  through  a  3-day  period,  changing  surface 
temperatures  with  iterations.  The  constant  temperature  is  a  valid  assumption 
since  at  the  starting  time,  the  temperature  is  nearly  constant  with  depth. 

The  finite  element  analysis  yielded  the  tautochrones  shown  in 
Figure  48  for  the  third  day.  The  convergence  is  quite  rapid  as  the  second  and 
third  day  tautochrones  are  nearly  the  same.  By  comparing  the  actual  tautochrones 
(Figure  46)  with  the  generated  values  (Figure  48),  close  agreement  -  within  3  or 
4°F  (1.8  to  2.2°C)  -  can  be  seen.  The  trends  in  the  curves  are  accurately 
portrayed  with  thermal  damping  effects  simulated  well. 

This  verification  has  only  considered  temperatures  above  freezing. 
The  program  does  not  accurately  account  for  the  water  phase  changes  associated 
with  freezing  and  thawing.  Provisions  can  be  made,  however,  to  deal  with  this 
complex  problem. 

3.3  VERIFICATION  FOR  MECHANICAL  LOAOS 

The  experimental  pavement  study  was  discussed  in  Chapter  II.  By  deter- 
mining the  geometry,  boundary  conditions,  and  material  properties  of  the  physi- 
cal model,  a  mathematical  model  can  be  gentrated  and  results  compared.  A 
comparison  of  anchored  pavement  with  the  conventional  slab  and  with  the  elevated 
anchored  slab  (representing  loss  of  contact  below  the  slab)  is  presented.  First, 
the  input  properties  of  different  materials,  used  in  analyses,  are  identified. 

The  properties  of  the  materials  utilized  in  the  experimental  study  are 
presented  in  Appendix  B.  Soil  properties  (E  and  v)  were  determined  by  several 
methods,  including  triaxial  compression,  constant  stress  consolidation,  and 
plate  load  test. 
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Figure  46.     Actual  recorded  tautochrones  at  Sta.  724+7.5  of  the 
Vormsund  test  road,  April  26,  1963  (Nordal ,  1973). 
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Table  1.     Material   properties  used  in  thermal   verification  program. 


MATERIAL 

THERMAL 
CONDUCTIVITY 

(BTU/HR-IN-F0) 

MASS  SPEC] 
HEAT 

FIC 

DENSITY 

(BTU/LBM- 

F°) 

(LB/IN3) 

Bitumen 

0.069 

0.245 

0.082 

Cr.  Stone 

0.093 

0.267 

0.075 

Gravel 

0.082 

0.226 

0.075 

Old  Gravel 
Pavement 

0.082 

0.226 

0.075 

Clay 


0.059 


0.345 


0.067 


1   BTU/HR-IN-°F  = 


0.099 


1.8°C  +  32 


CAL-KG 
HR-CM-°C 


1   BTU/LBM-°F  =  1t8gc5f32  ^T 


1   LB/IN3  =  27.7  —. 
CMV 


60 


TEMPERATURE  (8F) 


Ul 

o 


NUMBERS  INDICATE  HOUR  OF  THE  DAY 


1  in.  *  2.54  cm 
X   «  (°F  -  32)/1.8 


Figure  48.  Tautochrones  generated  by  the  finite  element  model 
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It  is  reasonable  to  assume  that  the  plate  load  test  will  provide  a 
better  in-situ  value  of  the  modulus  of  elasticity.  The  plate  load  test  was 
conducted  in  the  sampling  area.  This  area  did  not  undergo  many  preloadings  as 
the  main  experimental  area.  The  preloadings  will  increase  the  modulus.  The 
adopted  modulus  for  the  analytical  comparison  was  350  psi  (60  kg/cm  ) .  The  input 
values  are  summarized  in  Table  2. 

Eight  node  brick  elements  were  used  in  the  finite  element  analysis 
(3-0  analysis).  In  the  case  of  AS62  and  S62,  interface  elements  of  infinitesimal 
length  were  used  between  the  slab  and  the  subgrade  soil.  Thus  the  effect  of  the 
uplifting,  when  the  slab  looses  contact  with  the  subgrade  soil,  was  taken  into 
consideration.  Modeling  the  steel  container  was  considered  necessary  to  avoid 
great  errors,  as  found  by  preliminary  analysis. 

3.3.1  Results  of  Analysis 

The  finite  element  models  of  simply  supported  anchored  slab  (SAS62) 
and  simply  supported  conventional  slab  (SS62)  are  shown  in  Figure  49.  A  line 

load  of  2.95  psi  (0.21  kg/cm)  was  applied  1  in.  (2.54  cm)  offset  of  the 

transverse  centerline.  From  the  "beam  theory"  the  following  relation  can  be 
computed: 


ASAS62     :SAS62     n/[ 
ASS62     ISS62 

where 

A  denotes  deflection  under  load 

and 

I  denotes  moment  of  inertia 

The  equivalent  moment  of  intertia  for  the  anchored  slab  (SAS62)  is  5.28  in. 
(219.78  cm4)  and  for  the  conventional  slab  (SS62)  is  0.22  in.4  (9.16  cm  ).  The 
neutral  axis  of  the  cross  section  of  the  anchored  slab  is  0.647  in.  (1.65  cm) 
below  the  top  surface  of  the  slab. 

The  longitudinal  deflections  at  the  centerline  are  plotted  in 
Figure  50  showing  the  experimental  deflection,  the  deflection  obtained  by  finite 
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Table  2.  Material  properties  used  in  stress  verification  program. 

POISSON'S  RATIO 

0.35 
0.33 
0.40 


MATERIAL 

YOUNG'S  MODULUS 

(psi) 

Aluminum 

10  x  106 

Steel 

30  x  106 

Soil 

(SUty  Clay) 

850 

1  psi  =  0.0703 

2 

kg/  cm 

63 


a)  Finite  element  model  of  SAS62. 


SAS62  represents  a  simply 
supported  anchored 
slab 


b)Finite  element  model  of  SS62. 


SS62  represents  a  simply 
supported  convention- 
al slab 


Figure  49. 


Finite  element  models  of  SAS62  and  SS62  for  stress 


analysis  verification. 
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Figure  50.  Longitudinal  deflections  at  center  for  comparison  of 
experimental  models,  finite  element  models,  and  models  using 
the  beam  theory. 
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element  method,  and  that  obtained  by  beam  theory.  The  difference  between 
the  results  from  finite  element  method  and  the  experimental  values  is  about 
7%   for  the  simply  supported  anchored  slab  (SAS62)  and  about  4.5%  for  the 
simply  supported  conventional  slab  (SS62).  The  difference  between  the  beam 
theory  approach  and  the  experimental  value  is  about  Q%   for  the  anchored  slab 
and  about  8.5%  for  the  conventional  slab. 

Figure  51  shows  the  longitudinal  deflections  at  the  centerline 
of  the  simply  supported  anchored  slab  and  simply  supported  conventional  slab, 
using  the  finite  element  approach.  The  finite  element  models  for  anchored 
slab  on  subgrade  (AS62)  and  conventional  slab  on  subgrade  (S62),  are  shown 
in  Figure  52.  The  finite  element  model  for  an  elevated  slab  (ASE52)  is  shown 
in  Figure  53.  A  750-lb.  (340.2  kg)  load  was  applied  at  position  1  (Figure  8) 
in  the  experiment  sequence  AS62  and  S62;  a  200-lb.  (90.72  kg)  load  was  applied 
at  position  1  in  the  sequence  ASE62. 

Figure  54  shows  the  surface  deflection  contours  for  the  sequence 
S62  with  a  maximum  deflection  of  0.1  in.  (2.54  mm)  under  the  load.  More 
than  half  of  the  slab  looses  contact  from  the  subgrade  soil,  with  a  maximum 
uplift  of  0.012  in.  (0.30  mm  ).  The  vertical  stress  and  strain  contours 
are  given  in  Figure  55,  with  a  maximum  stress  of  9  psi  (0.63  kg/cm  )  and 
a  maximum  strain  of  0.005  in. /in.  (cm/cm)  underneath  the  load. 

Figure  56  shows  the  surface  deflection  contours  for  the  sequence 
AS62  with  a  maximum  deflection  of  0.031  in.  (0.8  mm)  in  the  right  anchor 
(closest  to  the  load).  No  serious  uplift  of  the  anchored  slab  is  noted. 
Only  the  left  and  right  edges  are  lifted,  with  a  maximum  uplift  of  0.004 
in.  (0.1  mm).  The  vertical  stress  and  strain  contours  in  the  subgrade  are 
given  in  Figure  57  with  a  maximum  stress  of  2  psi  (0.14  kg/cm  )  and  a  maximum 
strain  of  0.002  (in. /in.  or  cm/cm)  underneath  the  right  anchor.  The  stresses 

in  the  subgrade  immediately  under  the  slab  are  reduced  considerably  with 

2  2 

a  maximum  stress  of  1  psi  (0.07  kg/cm  ),  compared  with  9  psi  (0.63  kg/cm  ) 

in  conventional  slab. 

Figure  58  shows  the  surface  deflection  contours  for  the  sequence 
ASE62  with  a  maximum  deflection  of  0.011  in.  (0.28  mm)  at  the  right  edge 
of  the  centerline  (near  the  load),  with  no  serious  uplift.  Only  the  left 
far  edge  is  lifted,  with  a  maximum  uplift  of  0.001  in.  (0.03  mm).  The  stress 
and  strain  contours  are  plotted  in  Figure  59,  with  a  maximum  stress  of  1  psi 
(0.07  kg/cm  )  and  a  maximum  strain  of  0.0008  (in. /in.  or  cm/cm)  underneath 
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a)  Finite  element  model  of  AS62. 


AS62  represents  62  in.  long 
anchored  slab  on  soil 
subgrade 


b)  Finite  element  model  of  S62. 


S62  represents  62  in.  long 
slab  on  soil  subgrade 


Figure  52.  Finite  element  models  of  AS62  and  S62  for  stress 
analysis  verification. 
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ASE62  represents  62  in.  long 
anchored  slab  elevated 
0.5  in.  on  soil  subgrade 


Figure  53.  Finite  element  model  of  ASE62  for  stress  analysis 

verification. 
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1   in. ' -  2.54  cm 


S62  represents  62  in.   long 
slab  on  soil   subgrade 


Figure  54.     Surface  deflection  contours   (in.)  for  750-1 b.    (340kg) 

load  at  position  1   of  S62. 


70 


1  psi  =0.07  kg/cm 
a)  Vertical  stress  (psi)  contours. 


S62  represents  62  in.  long 
slab  on  soil  subgrade 


1  micro  in. /in.  -   0.001  cm/ cm 
b)  Vertical  strain  (microin./in.)  contours 


Figure  55.  Vertical  stress  and  strain  contours  for  750-lb 
(340  kg)  load  at  position  1  of  S62. 
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1   in.  =  2.54  cm 

AS62  represents  62  in.   long 
anchored  slab  on  soil 
subgrade 


Figure  56.     Surface  deflection  contours   (in.)  for  750-1 b. 
(340  kg)   load  at  position  1   of  AS62. 
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1   psi  -  0.07  kg/cm 


a)  Vertical  stress  (psi)  contours. 

AS62  represents  62  in.   long 
anchored  slab  on  soil 
subgrade 
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1  micro  in. /in.  ■  0.001  cm/cm 


b)  Vertical   strain   (micro  in. /in.)  contours 


Figure  57.     Vertical   stress  and  strain  contours  for  750-1 b 
(340  kg)   load  at  position  1   of  AS62. 
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1   in.'  »  2.54  cm 


ASE62  represents  62  in.  long 
anchored  slab  elevated 
0.5  in.  on  soil  subgrade 


Figure  58.     Surface  deflection  contours  (in.)  for  200-1  b 
(91   kg)  load  at  position  1  of  ASE62. 
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1   psi  a  0.07  kg/ cm 
a)  Vertical  stress  (ps1)  contours 


1  micro  1n./in.  ■  0.001  cm/cm 
b)  Vertical  strain  (micro  in. /in)  contours. 

ASE62  represents  62  in.  long 
anchored  slab  elevated 
0.5  in.  on  soil  subgrade 


Figure  59.     Vertical  stress  and  strain  contours  for  200-1 b 
C91   kg)  load  at  position  1  of  ASE62. 
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the  right  anchor  (near  the  load).  The  central  anchors  carry  23%  less  than 
the  right  anchors,  while  the  left  anchor  carries  49%  less  load  than  the  right 
anchor. 

The  finite  element  results  are  compared  with  the  experimental  analy- 
sis. For  sequence  S62,  the  transverse  surface  deflections  of  the  finite 
element  model  along  the  transverse  center  line  Tl-Rl  for  750-1  lb.  (342  kg) 
load  at  position  1,  are  plotted  in  Figure  60.  The  difference  is  about  20%. 
The  conventional  slab  was  loaded  after  the  anchored  and  elevated  anchored 
pavement  loading  was  completed.  It  is  apparent  that  many  loadings  and  unload- 
ings  have  strengthened  the  soil.  This  is  probably  one  of  the  major  factors 
contributing  to  this  difference  between  experimental  and  analytical  values. 
In  Figure  61,  the  transverse  and  longitudinal  centerline  surface  deflections 
of  the  finite  element  model  of  AS62  for  750-1  lb.  (342  kg)  load  are  plotted. 
The  difference  is  about  11%.  The  longitudinal  and  the  transverse  surface 
deflections  of  the  finite  element  model  of  ASE62  are  shown  in  Figure  62 
for  a  300-lb.  (157  kg)  load  at  position  1.  The  difference  of  the  transverse 
deflection  (Tl-Rl)  between  the  finite  element  and  the  experimental  results 
of  ASE62  is  about  18%.  In  general,  the  trends  of  observed  deflections  of 
the  model  tests  agree  well  with  the  analytical  (finite  element)  results. 

Although  the  results  of  only  a  few  of  the  many  experiments  conducted 
are  presented  in  this  report,  the  inferences  are  based  on  a  study  of  all 
the  experimental  results. 
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1  in, 
1  lb. 
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S62  represents  62  in.  long 
slab  on  soil  subgrade 


Figure  60.  Transverse  surface  deflection  along  line  Tl-Rl, 
load  position  1  of  S62. 
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1   in.   a  2.54  cm 
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a)  Transverse  surface  deflection  along  line  Tl-Rl 
load  position  1,  AS62. 
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b)  Longitudinal  surface  deflection  along  line  L2-N2,  load 
load  position  1,  AS62. 


Figure  61.     Transverse  and  longitudinal   surface  deflection  along 
lines  Tl-Rl   and  L2-N2,~load  position  1,  AS62. 
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b)  Longitudinal   surface  deflection  along  line  L2-N2, 
load  position  1 ,  ASE62. 

Figure  62.     Transverse  and  longitudinal   surface  deflection  alona 
lines  Tl-Rl   and  L2-N2,   load  position  1,  ASE62. 
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CHAPTER  IV 
ANALYSIS  OF  THE  PAVEMENT  SYSTEMS 

4.1  INTRODUCTION 

The  goal  of  the  analysis  was  to  compare  the  behavior  of  the  proposed 
pavement  system  (the  anchored  pavement)  with  that  of  a  conventional  pavement. 
The  physical  configuration  of  each  pavement  has  been  shown  in  Chapter  I 
(Figure  1).  The  results  of  the  analysis  are  divided  into  two  parts  - 
environmental  loading  effects  and  mechanical  loading  effects.  The  results  are 
later  superimposed.  A  pavement  joint  is  also  analysed  for  suitability  in  the 
proposed  anchored  pavement. 

4.2  MATERIAL  CHARACTERIZATION 

4.2.1  General  Behavior  of  Pavement  Systems 

Pavement  systems  generally  do  not  exhibit  the  same  behavior  throughout 
the  year.  In  areas  prone  to  freezing  or  freeze/thaw  cycles,  pavement  behavior  is 
often  radically  different  under  frozen  and  unfrozen  conditions,  especially  in 
frost-susceptible  soils  such  as  clays  and  silts.  Carpenter  (1976)  and  others 
have  attempted  to  isolate  gross  climatic  conditions  in  the  United  States  using 
descriptors  such  as  Wet-Freeze/Thaw,  Dry-No  Freeze  (Figure  P3).  By  rationally 
determining^  the  extreme  conditions  which  occur  with  a  high  enough  probability,  it 
is  possible  to  design  a  pavement  system  to  withstand  these  extremes  as  well  as 
more  normal  conditions. 

As  an  example,  consider  the  results  reported  by  Moore  et  al.  (1969)  for 
surface  deflections  of  a  flexible  pavement  in  a  wet,  seasonal  freeze-thaw  area. 
Figure  64  presents  a  trend  in  surface  deflection  characteristics  throughout  a 
typical  year.  Four  primary  periods  exist  in  the  yearly  cycle:  1)  period  of  deep 
frost  characterized  by  low  deflections  due  to  increased  rigidity  of  the  pavement 
as  well  as  increased  bearing  capacity  and  modulus  of  subgrad.-*  reaction  of  the 
soil;  2)  period  of  rapid  strength  loss  during  the  spring  thaw,  which  gives  rise 
to  a  corresponding  rapid  increase  in  surface  deflections;  3)  period  of  rapid 
recovery  in  late  spring;  and  4)  period  of  gradual  slow  strength  recovery  leading 
back  to  the  frost  condition.  There  would  probably  be  a  rather  rapid  strength 
recovery  after  the  first  permanent  frost  preceding  the  period  of  low  deflections, 
thus  starting  the  cycle  again.  Shapes  of  deflection  basins  at  certain  key  times 
can  be  seen  in  Figure  65.  The  almost  uniform  deflection  in  the  frost  period  and 
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DISTANCE   FROM  LOAD 
1  2 


(FT) 


O  MID-WINTER 

O  MID-FALL 

a  EARLY   SPRING 


1   ft.    =   30.48    an 
1  mil  ■   0.001   in. 


Figure  65.     Comparison  of  deflection  basins  for  various  times  of 
year  under  1000  pound  load  (453.6  kg)  of  dynaflect 
test  apparatus  (Moore  et  al .   1969). 
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the  more  localized  deformation  during  the  period  of  subgrade  weakening  in  early 
spring  can  be  easily  noted. 

One  obvious  critical  period  in  the  structures'  design  life  is  the 
cyclic  appearance  of  these  periods  of  rapid  strength  loss.  This  problem  is 
aggravated  by  the  upward  flow  of  moisture  in  both  liquid  and  vapor  form  as  soil 
freezes  due  to  flow  along  thermal  gradients  and  also  due  to  subpressures  at  the 
freezing  front.  When  the  thaw  period  arrives,  a  greater  than  normal  water 
content  exists  in  the  upper  layer  of  the  subgrade,  and  the  soil  is  greatly 
weakened.  Figure  66  shows  the  typical  seasonal  change  in  moisture  content  near 
the  surface  of  various  classes  of  soil.  In  general,  the  water  content  increases 
until  the  thaw  period  then  decreases  gradually.  The  finer  grained  soils  (clays), 
because  of  lower  permeability,  exhibit  smaller  changes  in  moisture  throughout 
the  year.  Also,  during  summer,  as  temperatures  increase,  evapo-transpiration 
tends  to  decrease  moisture  content  until  the  freezing  period  when  the  cycle  is 
reinitiated. 

4.2.2  Material  Properties 

Two  basic  types  of  analysis  were  performed:  heat  transfer  analysis  and 
mechanical  stress  analysis.  The  finite  element  formulation  used  for  analysis 
requires  knowledge  of  certain  material  properties.  ANSYS  is  capable  of  handling 
temperature-dependent,  nonlinear  material  mechanical  behavior  and  temperature- 
dependent  thermal  properties  such  as  thermal  conductivity  and  heat  capacity. 

Due  to  the  extreme  variability  in  soil  properties  from  soil  to  soil  and 
with  variation  in  moisture,  temperature,  confining  pressure,  etc.,  it  was 
decided  to  choose  typical  extreme  values  for  properties  to  represent  limiting 
(best  or  worst)  cases  rather  than  attempt  a  parametric  study  of  pavement  response 
under  several  types  of  material  constituents. 

4.2.2.1  Thermal .  Heat  transfer  through  soils  is  primarily  a  conduction  phe- 
nomenon; however,  convection  and  radiation  do  occur.  Solar  radiation  and  con- 
vection through  air  and  liquids  at  the  soil  surface  as  well  as  convection  of  heat 
through  vapor,  soil,  and  liquid  interfaces  within  the  pore  structure  are 
possible.  In  the  basic  three-dimensional  heat  transfer  equation  for  conduction, 

ell  =  -Ifk  II)  +  _l(k  ^1)  +  -i(k  H )  4 

u3t   3x^  x3x;  dyK  ydyJ       ZzK  zdzJ 
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two  material  properties  are  required:  heat  capacity,  c,  and  conductivity,  k. 
Heat  capacity  is  defined  as  the  quantity  of  heat  necessary  to  produce  a  'unit 
change  in  temperature  for  a  unit  mass.  Because  of  the  multi -phase  character  of 
soil  (vapor,  liquid,  minerals),  the  total  heat  capacity  is  usually  a  weighted  sum 
of  the  components  such  that 

c  =  Z  gi  ci  5 

4.L. 

where  g.   is  the  mass  of  the  i —  component  and  c,   is   its  heat  capacity  per  unit 
mass.     A  specific  relation  can  be  made  for  frozen  and  unfrozen  soils: 

WT 
:v  =  Vms  +  IDff  ^danw  6 


T>0°C:    C   =  V,C-  +  T*7 


o  WT        WIT 

T<0  C:   cy  =  Ydcms  +  m   Ydcmw  +  m  Y(jcmi 

7 

,  /Aw/ 100 A   . 
+  (-^V  dL 

where  y.   is  soil  dry  density;  c  is  volmetric  heat  capacity;  L  is  latent  heat  of 
fusion  of  water;  c_s,  c  ,  and  c  •  are  mass  heat  capacities  of  soil,  water,  and 
ice,  respectively;  Wj  and  w,j  are  moisture  contents  of  unfrozen  and  frozen  water 
at  temperature  T,  respectively;  and  the  term  in  parenthesis  is  the  gradient  of 
unfrozen  water  content  versus  temperature  at  temperature  T  (Williams,  1973).  The 
heat  capacity  of  moist  soil  as  it  freezes  is  thus  augmented  by  the  latent  heat  of 
fusion  of  water,  which  produces  a  large  spike  in  the  c  versus  T  relation,  often 
20  to  200  times  as  great  as  unfrozen  or  well-frozen  soil  heat  capacities,  depend- 
ing upon  water  content.  This  spike  generally  appears  over  a  small  range  of 
temperatures,  usually  less  than  1°F  (0.56°C)  (Figure  67).  Some  typical 
volumetric  heat  capacities  for  soils  are  presented  in  Figure  68  showing  varia- 
tion in  c  with  temperature  as  a  soil  freezes  and  as  it  thaws.  Table  3  presents 
values  for  soils  and  structural  materials  that  are  used  in  the  finite  element 
models.  Notice  that  asphaltic  concrete  and  plain  concrete  can  have  peaks  if  some 
moisture  is  present  in  pores. 

Conductivity  is  a  measure  of  the  ease  with  which  heat  can  flow  through 
a  material.  Some  typical  variations  of  thermal  conductivity  with  temperature  for 
soils  and  structural  materials  are  shown  in  Figure  69.  Kersten  (1949)  has  made 
extensive  laboratory  tests  to  measure  conductivity  of  soils  and  has  developed 
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Figure  67.     Conceptual  dependence  of  thermal   properties  on 
temperatures   (Frivik  et  al.,  1977;  Comini   et  al . ,   1974). 
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Figure  68.     Volumetric  heat  capacities  of  various  soils  near 
freezing  temperatures  (Williams,  1973). 
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Figure  59.     Variation  of  thermal   conductivities  with  tamoerature  (Friv'ik 
at  al.,   1977;  Seville,   1963;  Sucac,   1975;  Hoe.kstra,   1969). 
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Table  3.  Typical  heat  capacities  of  road  materials  (BTU/lb.-F°) 

(Frivick  et.  al.,  1977). 


TEMPERATURE  (°F) 

MATERIAL 

-22.00 

30.38 

31.98 

32.00 

122.00 

Bitumen 

0.22 

0.22 

10.1 

0.26 

0.29 

Plain  Concrete 

0.16 

Coarse  Gravel 

0.18 

0.18 

5.06 

0.19 

0.22 

Sandy  Gravel 

0.23 

0.23 

29.0 

0.31 

0.33 

Clay 

0.22 

0.22 

33.4 

0.30 

0.32 

X  -  (°F  -  32)/1.8 
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empirical  relations  for  thermal  conductivity  of  frozen  and  unfrozen  sands  and 
clays  based  on  water  content  (Figure  70).  His  general  observations  are  as 
follows: 

1.  k  of  unfrozen  soil  increases  slightly  with  temperature. 

2.  At  optimum  moisture  content  for  compaction,  the  k  below  freezing 
averages  about  17%  greater  than  above  freezing. 

3.  At  a  constant  moisture  content,  as  density  increases,  k  increases. 

4.  At  constant  density,  as  moisture  increases,  k  increases  up  to 
saturation. 

5.  For  saturated  unfrozen  soils,  k  decreases  as  density  decreases. 
For  frozen  soils,  there  is  no  clear  relation. 

6.  k  varies  with  texture.  For  a  given  density  and  water  content,  k  is 
greater  in  coarser  soils  (gravel  and  sand)  than  in  fine  soils 
(silty  loam,  clay). 

7.  k  varies  with  mineral  content.  Quartz  has  a  high  k,  kaolin  a  low 
value,  and  feldspars  an  intermediate  value. 

Thermal  conductivity  has  also  been  linked  to  strength,  as  shown  in  Figure  71. 
This  relation  is  rather  crude  and  is  a  composite  of  several  curves  for  individual 
materials  (clays,  sands,  etc.). 

Convection  at  the  surface  has  also  been  considered.  Ouffie  and  Beckman 
(1974)  report  a  relation  for  the  film  coefficient  of  forced  convection  heat  loss 
from  flat  plates  exposed  to  wind: 

h,  -   5.7  +  3.8v  8 

w 

?   r» 
where  h  is  the  film  coefficient  or  heat  transfer  coefficient  in  W/m  -  C  and  v 
w 

is  wind  velocity  in  m/sec. 

4.2.2.2.  Mechanical.  Mechanical  strength  properties  of  materials  within  pave- 
ment systems  also  vary  with  temperature.  This  temperature  dependence  is 
important  during  a  thermal  stress  analysis  and  reflects  material  behavior  of  the 
system  that  may  be  greatly  different  over  different  temperature  ranges.  The 
temperature-dependent  strength  properties  used  in  this  investigation  were 
elastic  modulus  (for  a  linear  analysis),  Poisson's  ratio,  and  the  coefficient  of 
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Figure  70.  Variation  in  thermal  conductivity  of  soils  with 
moisture  content,  (for  100  PCF  dry  density) 
(Kerston,  1949). 
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Figure  71.  General  correlation  between  modulus  and  thermal 
conductivity  in  frozen  soils  (Jumikis,  1966). 
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linear  expansion  of  the  material  under  a  uniform  temperature  gradient.   In  a 
nonlinear  inelastic  analysis,  strength  is  input  as  a  series  of  piecewise  linear 
stress-strain  curves  with  some  initial  elastic  modulus. 

Because  of  the  complex  multi-material  and  multi-phase  constituency  of 
soil,  it  is  difficult  to  characterize  behavior  over  a  complete  temperature  regime 
with  only  a  modulus  and  Poisson's  ratio.  As  moist  soil  freezes,  all  water  within 
the  soil  system  does  not  freeze  at  once.  Jumikis  (1966)  has  described  the 
freezing  process  in  soils  by  first  considering  the  types  of  water  within  soil: 
free   water,  which  exists  in  pores  in  bulk  and  is  not  bound  to  the  soil  particles; 
capillary  water,  which  is  held  in  sufficiently  small  pore  structures  by  capillary 
action;  film  water  surrounding  particles  in  various  layers;  and  hygroscopic 
moisture  actually  bound  to  the  soil  (film  water  and  hygroscopic  water  are 
adsorbed  water).  These  forms  of  moisture  usually  exist  simultaneously  and  have 
different  freezing  points,  resulting  in  quite  complex  microscopic  behavior. 
Free  water  freezes  first  at  32°F  (0°C).  Capillary  water  freezes  at  about  29.5°F 
(-1.4  C)  and  freezing  precedes- through  the  less  stressed  film  layers.  Film  water 
is  stressed  by  the  electrostatic  attraction  between  the  water  di poles  and  the 
charged  surface  at  the  soil  particle,  so  film  water  close  to  the  particle  surface 
is  more  highly  stressed  and  has  a  lower  freezing  point. 

A  simplified  view  of  this  action  is  depicted  in  Figure  72.  An 
arbitrary  division  of  the  film  into  four  layers  illustrates  the  freezing  process. 
Layer  1  freezes  first,  then  layer  2,  etc.  As  the  stressed  water  freezes,  the 
bonding  forces  between  water  molecules  are  broken,  and  they  begin  to  leave  the 
film  water  and  join  growing  ice  crystals.  Because  free  water  freezes  first,  film 
water  tends  to  migrate  to  the  freezing  bulk  water  ice  crystals  at  and  behind  the 
freezing  front. 

To  simulate  the  mechanical  behavior  due  to  this  gradual  accumulation 
and  growth  of  ice  in  the  pore  structure  of  moist  soils,  a  modulus  versus  tempe- 
rature relation  was  developed.  Figure  73  depicts  possible  soil  behavior  as  film 
water  freezes,  binding  soil  particles  in  a  matrix  of  frozen  water.  The  relation 
is  conceptual  and  would  vary  greatly  from  soil  to  soil.  Coarse-grained  soil,  in 
which  a  higher  percentage  of  water  is  free  water,  would  tend  to  have  a  more  rapid 
strength  increase  as  the  temperature  dropped  below  freezing  due  to  a  smaller 
unfrozen  water  content  (Figure  74).  Finer  grained  soils,  such  as  clays,  with 
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Figure  72.  Conceptual  diagram  illustrating  the 
freezing  of  soil  moisture. 
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Figure  73.     Effect  of  temperature  on  modulus  of  moist,  fine-drained  soil 
due  to  "cementation"  of  particles  by  ice  (Tystovich,  1975). 
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high  specific  surface  areas,  would  probably  have  the  gradual  strength  increase 
shown  in  Figure  73  (Tsytovich,  1975). 

The  modulus  of  elasticity  of  ice  varies  with  temperature  and  type  of 
ice,  but  is  generally  in  the  range  of  300,000  to  1,500,000  psi  (21,000  to  105,000 
kg/on  ) .  With  this  in  mind,  a  crude  relation  between  modulus  and  temperature  may 
be  drawn  based  upon  a  weight  relation  such  that 

£T  =  (wyj/lOOJEj  +  (l-wTI)/100)  Es  9 

where  Ep  Et,  and  E  are  moduli  of  the  total  soil  system,  ice,  and  unfrozen  soil, 
respectively,  and  Wj,  is  the  percentage  of  ice  in  the  system.  This  relation, 
however,  does  not  consider  the  matrix-particle  effect  attributed  to  the  altera- 
tion of  the  gross  structure  of  the  unfrozen  soil /water  system. 

Static  soil  strength  for  temperatures  above  freezing  is  not  so  vari- 
able. Typical  compressive  strengths  for  clays  at  different  temperatures  and 
water  contents  are  shown  in  Figure  75.  For  a  given  moisture  content,  clays 
decrease  in  strength  with  increasing  temperature  above  70 °F  (21°C),  decreasing 
approximately  30%  in  strength  for  a  change  of  75°F  (42°C).  Modulus  is  however 
less  variable  (Figure  76).  Poisson's  ratio  for  ice  is  about  0.33  to  0.38; 
similar  ranges  exist  in  unfrozen  soil.  Thus,  it  was  assumed  that  a  near  constant 
value  of  between  0.35  and  0.45  for  Poisson's  ratio  could  be  used  over  an  expected 
range  of  temperatures. 

The  coefficient  of  cubic  expansion  for  water  as  a  function  of  tem- 
perature is  shown  in  Figure  77.  Because  it  is  a  solid,  ice  shows  relatively 
little  volumetric  expansion  under  temperature  gradients.  Water  tends  to  rise 
linearly  from  the  freezing  point  until  near  the  vapor  phase. 

Structural  materials,  such  as  asphaltic  concrete  (AC),  Portland  cement 
concrete  (PCC),  and  steel  fiber  reinforced  concrete  (SFRC),  also  have  tempe- 
rature-dependent material  properties,  although  under  different  mechanisms.  AC 
is  a  viscoelastic  material  that  flows.  As  the  temperature  increases,  its  vis- 
cosity decreases  and  flow  increases,  leading  to  a  structural  weakening.  Figure 
78  depicts  this  phenomenon  by  a  rapid  decrease  in  modulus  at  temperatures  above 
50  F  (10°C).  Poisson's  ratio  for  AC  generally  increases  with  temperature  and  is 
limited  to  a  value  of  0.5,  that  of  an  incompressible  fluid. 
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Figure  76.     Effect  of  temperature  on  the  stress-strain  behavior  of 

a  clay  (Murayama,  1969). 
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Figure  78.  Strength  properties  of  asphaltic  concrete  (AC)  as  a 
function  of  temperature  (Nair  and  Chang,  1973). 
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Compressive  strength  of  PCC  does  not  vary  greatly  at  ambient  tem- 
peratures above  about  50°F  (10°C)  (Figure  79),  but  does  increase  at  freezing 
temperatures,  probably  through  a  mechanism  of  moisture  freezing  within  pores  in 
the  concrete-  Elastic  modulus  (E  )  for  PCC  has  been  extrapolated  using  the  ACI 
formula: 

Ec  =  57,500  /f'c  10 

to  provide  the  temperature-dependent  relation  shown  in  Figure  80.  Properties  of 
SFRC  were  then  extrapolated  using  the  same  relative  modulus  curve  and  a  70  F 
(21°C)  modulus  value  of  9  million  psi  (633,000  kg/cm2). 

The  material  properties  used  in  the  finite  element  modeling  procedure 
are  summarized  in  Table  4.  Linear  interpolation  and  extrapolation  for  inter- 
mediate values  are  automatically  used  in  the  computer  program. 

4.3  ENVIRONMENTAL  LOADS 

4.3.1  General 

Environmental  loadings  on  pavement  systems  are  defined  in  this  study 
as  consisting  of  thermal  loads  generated  by  radiation  and  convection  of  heat  at 
the  surface,  freeze/thaw  action,  localized  ice  lens  formations,  moisture  migra- 
tion, salt  intrusion,  and  factors  that  cannot  be  directly  solved  by  the  finite 
element  method,  such  as  the  effects  of  vegetation  and  chemical  attack.  Heat 
transfer  and  moisture  migration  can  be  analyzed  directly  by  finite  element 
analysis.  Forms  of  chemical  attack,  the  effect  of  changes  in  moisture  dis- 
tribution, etc.,  can  be  analyzed  indirectly  by  considering  the  effect  of  these 
actions  on  the  properties  of  the  system,  then  performing  a  stress  or  thermal 
analysis. 

4.3.2  Characterization  of  Loads 

The  primary  thrust  of  the.  investigation  centers  on  the  effects  of 
thermal  loads.  To  fully  understand  the  nature  of  thermal  inputs  into  a  pavement 
system,  a  knowledge  of  the  heat  balance  at  the  pavement  (or  soil)  surface  is 
required.  For  an  infinitesimal  thickness  of  surface  it  is  possible  to  express 
the  balance  of  heat  flow,  both  into  and  out  of  the  surface.  Figure  81  depicts 
the  balance  at  the  surface,  including  heat  generated  within  the  soil  and  heat 
input  at  the  surface  as  convection  and  radiation  loads.  External  sources  of  heat 
are  subdivided  into  four  categories.  First  is  radiation,  perhaps  the  largest 
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Table  4.    Material  properties  used  in  the  finite  element  model. 
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Figure  81.     Heat  balance  at  the  pavement  surface  (Noss,  1973). 
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single  heat  source  and  most  important  factor  in  the  heat  balance  at  the  surface. 
The  short-wave  radiation  source,  the  sun,  contributes  direct  and  diffuse  radia- 
tion (Figure  82).  Long-wave  radiation  from  clouds  and  water  vapor  and  long-wave 
radiation  directed  outward  from  the  surface  complete  the  sources. 

Secondly,  direct  heat  flow  due  to  convection  and  conduction  at  the 
surface  both  add  and  withdraw  heat  from  the  surface.  Convection  in  pavement 
systems  consists  almost  entirely  of  forced  convection,  heat  transferred  through 
a  film  of  air  moving  across  the  surface.  Equation  3  gives  a  simple  relation  to 
calculate  the  coefficient  of  heat  transfer  or  film  coefficient  for  a  flat  hor- 
izontal surface  in  wind.  Convection  can  be  modeled  with  convection  elements 
available  in  the  program  that  requires  as  input  the  film  coefficient  and  con- 
tributing area.  This  convection  link  is  a  linear  element  with  two  nodes.  One 
node  receives  the  temperature  input,  and  the  other  is  connected  to  the  surface. 
Conduction  is  also  possible  through  contact  with  materials  of  different  tem- 
perature at  the  surface. 

Phase  changes  constitute  another  important  source  of  heat  at  the 
surface.  As  water  freezes,  it  releases  a  quantity  of  heat  equal  to  the  latent 
heat  of  fusion.  This  heat  is  dissipated  at  the  surface  as  surface  ice  and  snow 
freezes.  Conversely,  as  surface  ice  thaws,  it  requires  a  heat  equal  to  the 
latent  heat  of  fusion  to  convert  from  solid  to  liquid.  The  latent  heat  of  fusion 
of  water  at  atmospheric  pressure  is  80  cal/gm. 

Evaporation  and  condensation  are  also  phase  change  processes,  but  from 
liquid  to  vapor.  Here  the  important  quantity  is  the  latent  heat  of  vaporization, 
which  for  water  at  atmospheric  pressure  is  540  cal/gm. 

Finally,  heat  exchange  on  the  surface  between  precipitation  and  the 
surface  is  possible.  Both  warm  (rain)  and  cold  (snow,  sleet)  precipitation  exist 
(depending  upon  locality),  and  heat  exchange  occurs  if  there  is  a  temperature 
differential  between  the  precipitation  and  the  surface. 

Internal  heat  sources  include  geothermic  heat,  from  deep  within  the 
ground,  latent  heat  attributed  to  phase  changes  within  the  soil  moisture,  and 
volumetric  heat  associated  with  the  temperature  of  the  system,  in  particular  any 
gradient  that  exists  between  the  soil  temperature  some  distance  from  the  surface 
and  close  to  the  surface. 
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The  surface  heat  balance  is  an  equilibrium  state;  however,  the  quan- 
tities are  seldom  constant.  Surface  temperature  is  not  constant  either;  it  is 
closely  linked  to  ambient  air  temperature  and  follows  the  same  trends.  As  air 
temperature  increases,  surface  temperature  increases  with  very  little,  if  any, 
phase  difference  or  time  lag.  In  winter  in  regions  with  frost  conditions,  there 
is  very  little  difference  between  air  and  surface  temperature  (Figure  83);  how- 
ever, during  summer  there  can  be  a  substantial  difference  -between  surface  and  air 
temperature,  often  as  much  as  10  to  20°F  (5  to  10°C),  as  indicated  in  Figure  34. 
(Surface  temperature  is  greater  during  the  day  due  to  radiation  effects). 

Surface  temperature  of  the  pavement  and  natural  ground  depends  on  many 
factors,  as  indicated  in  the  heat  balance  discussion.  The  temperature  due  to 
radiation  is  also  color  dependent;  pavements  with  dark  color  absorb  light  energy 
and  are  therefore  warmer.  This  condition  is  controlled  by  a  property  of  the 
surface  known  as  emissivity.  Ground  cover  conditions  can  greatly  influence 
surface  temperatures.  Vegetation  creates  a  rough  surface  so  that  convection  is 
different  from  that  existing  across  a  smooth  pavement  surface.  Snow  cover 
creates  a  reflective  surface  so  that  short-wave  radiation  is  reflected. 

The  temperatures  of  lower  layers  within  the  system  depend  not  only  upon 
the  heat  balance  at  the  surface,  but  also  on  heat  flow  to  and  from  lower  layers 
due  to  phase  changes,  conduction  of  heat  through  pore  fluid  and  intergranular 
contacts,  and  any  moisture  migration  along  thermal  gradients.  Typical  mean 
annual  changes  in  temperatures  at  various  depths  are  as  follows:  surface,  about 
50°F  (28°C);  at  12-in.  (30.5  cm)  depth,  about  10°F  (5°C);  at  24-in.  (61  cm) 
depth,  about  2°F  (1°C)  (after  Nordal,  1973).  The  soil  profile  at  depth  thus 
maintains  a  fairly  constant  temperature  while  surface  layers,  which  are  closer  to 
the  more  widely  variable  ambient  conditions,  possess  a  higher  degree  of 
variabil ity. 

A  typical  recorded  daily  change  of  temperature  in  the  desert  is 
presented  in  Figure  85.  The  surface  temperature  varies  to  a  large  degree,  being 
hottest  during  the  day  and  coldest  at  night.  Air  temperature  reflects  this 
change,  but  is  not  as  variable  and  does  not  reach  the  same  magnitude  of  tem- 
peratures as  the  surface.  Temperature  at  approximately  a  depth  of  24  in.  (61  cm) 
is  more  or  less  constant. 

Based  on  several  sets  of  similar  temperature  variations,  summer  and 
winter  daily  temperature  histories  were  simulated  by  the  models  shown  in 
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Figure  83.  Correlation  between  surface  and  air  temperatures 
(Venkatasubramanian,  1964). 
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Figure  36  and  Figure  87,  respectively.  The' summer  model  1s  rather  extreme  (upper 
boundary),  with  a  60°F  (33°C)  differential  between  maximum  and  minimum  surface 
temperatures.  Soil  surface  temperature  was  modeled  slightly  lower  than  the 
pavement  temperature,  but  this  1s  not  a  general  rule,  only  a  particular  example. 
Actual  temperatures  depend  on  ground  cover.  The  summer  condition  was  cycled  for 
2  days  and  rapidly  approached  a  steady  state  condition  (temperature  distribution 
was  nearly  the  same  for  corresponding  times  on  different  days). 

The  winter  model  attempts  to  generate  a  freeze/thaw  cycle  within  the 
system.  Most  of  the  history  is  of  subfreezing  temperatures.  This  model  was 
cycled  for  6  days,  but  convergence  to  a  steady  state  condition  as  defined  above 
and  a  constant  or  progressing  frost  line  1s  much  slower,  probably  because  of  the 
problems  associated  with  the  phase  change  of  water.  This  problem  and  a  solution 
are  discussed  later  1n  the  report. 

4.3.3  Description  of  Models 

The  first  analysis  performed  was  the  two-dimensional  heat  transfer. 
Two  modes  of  thermal  loading  were  considered:  simple  input  of  surface  tem- 
perature histories  (Figure  88)  and  a  surface  convection  analysis  (Figure  89)  in 
which  temperatures  were  Input  indirectly  through  convection  elements  to  simulate 
heating  and  cooling  of  the  surface  by  an  adjacent  film  of  air.  There  was  little 
difference  1n  the  resulting  temperature  distributions,  and  the  latter  model  was 
abandoned.  Because  Input  temperatures  are  derived  from  actual  measured  surface 
temperatures,  the  first  model  seems  justified. 

The  mesh  used  for  this  analysis  1s  shown  in  Figure  90.  Isoparametric 
elements  are  used  to  decrease  the  required  fineness  of  the  mesh  (number  of 
degrees  of  freedom)  and  still  maintain  accuracy.  Each  element  contains  a  three- 
by-three  Gaussian  grid  of  integration  points  for  the  Interpolation  of  tempera- 
tures (or  stresses)  within  the  element  from  the  calculated  nodal  temperatures 
(forces).  Each  node  has  one  degree  of  freedom  (temperature)  in  the  heat  transfer 
analysis,  and  two  degrees  of  freedom  per  node  (translation  in  the  plane)  in  the 
stress  analysis.  A  detailed  look  at  the  slab-soil  interface  (Figure  91)  reveals 
conduction  elements  (sliding  interface  elements  1n  the  stress  analysis)  to 
thermally  link  the  slab  and  soil.  To  aid  in  the  continuity  of  the  slab  and 
anchors,  a  temperature  degree  of  freedom  constraint  is  made  to  ensure  that 
corresponding  nodes  have  the  same  temperature.  This  modeling  technique  can  also 
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Figure  91.     Thermal  modeling  details 
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simulate  a  link  (dowels)  between  the  slab  and  anchors,  if  the  slab  and  anchors 
are  poured  separately. 

Boundary  conditions  used  in  the  model  to  simulate  an  infinite  depth  of 
the  soil  are  as  follows.  In  Figure  91  the  lower  boundary  is  fixed  at  a  constant 
temperature  (in  summer  50°F  (10°C)  and  in  winter  45°F  (7.2°C)).  This  boundary  in 
this  particular  model  is  about  50  ft.  (15.3  m)  below  the  surface  and  should 
accurately  simulate  no  heat  flow.  Lateral  boundaries  are  simulated  by  con- 
straining nodes  along  peripheral  nodal  lines  so  that  the  temperature  variation 
with  depth  is  the  same. 

4.3.4  Heat  Transfer  Analysis 

The  results  of  the  heat  transfer  analysis  used  to  generate  temperature 
distributions  within  the  pavement  and  soil  for  further  stress  analysis  are  pre- 
sented in  two  forms:  temperature  profiles  at  various  positions  and  times 
(tautochrones)  and  temperature  contours.  Typical  summer  and  winter  tautochrones 
for  the  anchored  and  conventional  pavements  are  shown  in  Figures  92  and  93, 
respectively.  Slight  differences  exist  in  the  profiles  under  the  anchors  due  to 
differences  in  conductivity  and  heat  capacity  between  soil  and  concrete.  The 
results  correlate  well  with  typical  field  data  (Figure  94)  in  trends  and  relative 
magnitudes  of  temperatures.  The  lag  in  temperature  peaks  with  increasing  depth 
was  accurately  portrayed  in  the  analysis  by  comparing  the  1800-hr.  and  2400-hr. 
tautochrones  with  the  1200-hr.  profile;  because  of  thermal  damping,  the  higher 
temperature  does  not  occur  at  lower  depths  until  much  later. 

Typical  temperature  contours  for  both  pavements  (anchored  and  con- 
ventional) have  been  studied  from  analyses  results,  but  they  are  not  included  in 
this  report.  The  temperature  gradients  are  not  as  steep  in  the  anchors,  due  to 
differences  in  material  properties  between  soil  and  concrete. 

4.3.5  Thermal  Stress  Analysis 

The  temperatures  generated  in  the  heat  transfer  analysis  were  input  in 
a  stress  model  and  the  resulting  stresses  and  strains  computed.  The  stress  model 
uses  sliding  interface  elements  that  cannot  sustain  tensile  loads  and  at  a 
certain  shear  stress  begin  to  slide  with  an  irrecoverable  displacement  as  shown 
in  Figure  95.  Actual  values  for  u,  the  coefficient  of  static  sliding  friction, 
and  k,  the  modulus  of  subgrade  reaction,  were  obtained  from  results  of  direct 
shear  tests  of  concrete  on  soil  reported  in  the  literature.  Relations  such  as 


121 


o  o  o  o 

©  ©  ©  © 

tf  (N  ffl  V 

©  cH  i-4  <S 


o 
u 


0) 


I 

> 


</1 

s 

s 

u 

2 

s 


i 

wo 


eg 


s 


cui)  aaYHO  moths  HXdaa 


122 


s-   i-     i.  i- 


O  O  O  O 

o  ©  ©  o 

VO   <M  00   T 

OH  fIN 

O    O  <3    0» 


o 

s 
u 

§ 


\  s 

— »  0 
CM 

ro  ^« 

in 

l     • 

CM 

b 
o     II 


e    iH 


1 

H 

EH 

Z 
> 

§ 


<u 


C 

<u 
u 


> 
ItJ 

« 
l/» 

c 
o 

i- 

o 


u 
<u 

■*-» 
c 

3 


en 


in 

CM 


o 
m 


m 


o 
o 


fUty    aaYHD  M013S   HidSQ 


123 


TEMPERATURE  (°F) 


a 

$ 

B 
O 

03 

SB 

04 

w 

Q 


0 
25 

50 

75        100 

12! 

^#^™ 

50 

p                        o  0600  hr. 
1          a  1400  hr. 

a  1900  hr. 

v  2300  hr. 

75 

• 

°C  =  (°F  -  32J/1.8 
1  in.  -  2.54  cm 

100 

Figure  94.  Tautochrones  in  natural  slope  on  a  summer 
day  (Goetez  and  Muller,  1969). 


124 


1 

FORC 

33, N 

/ 

T. 

/k 

(NO   TENSION) 

DISPLACEMENT 


a)  Normal   force/displacement  model. 


FORCE , T 
yT" 


T-. 


yT-- 


DISPLACEMENT 


b)  Shear  force/displacement  model. 
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those  found  in  Figure  96,  Teller  and. Sutherland  (1935),  and  Timms  (1964)  are 
typical.  In  most  instances,  k  depends  on  confining  pressure  (Figure  97)  or 
normal  stress,  a,  but  lower  bound  values  of  1000  lb. /in.  (178  kg/cm)  for  summer 
and  2000  lb. /in.  (356  kg/cm)  in  winter  were  used,  more  closely  resembling  clay 
behavior  with  increased  adhesion.  With  the  given  density  of  nodes,  this  corres- 
ponds to  a  modulus  of  subgrade  reaction  of  about  100  and  200  psi  (2.7  and  5.4 
kg/an  ),  respectively. 

With  this  type  of  interface  behavior,  it  is  possible  to  simulate  curl- 
ing of  the  slab  with  the  associated  loss  of  subgrade  contact  and  also  simulate 
penetration  of  the  slab  into  the  soil.  If  no  penetration  is  desired,  a 
constraint  limiting  the  compressive  deflection  of  the  slab  can  be  established  by 
choosing  a  sufficiently  large  value  for  k.  Surface  deflections  for  summer 
curling  conditions,  day  and  night  (Figure  98)  and  winter  conditions  (Figure  99) 
are  presented.  All  thermal  stress  calculations  are  based  on  a  zero  stress 
reference  temperature  of  68°F  (20°C),  the  assumed  temperature  of  concrete  when 
placed.  Daytime  deflections  are  concave  down,  due  to  the  warmer  surfaces; 
deflections  during  the  night  are  concave  up  because  of  the  cooler  surface. 
Deflections  are  somewhat  lower  in  winter  than  in  the  summer,  because  there  is  not 
as  great  a  temperature  differential  between  the  slab  top  and  bottom  surfaces. 
The  anchors  reduced  curling  near  the  pavement  edges  with  only  slight  curling  of 
the  shoulders.  Stresses,  however,  are  increased,  and  precautions  must  be  made  to 
increase  the  size  of  the  slab-anchor  joint  by  adding  fillets  to  lower  any  stress 
concentrations  that  may  occur. 

From  a  study  of  various  computer-generated  contour  plots  of  stress  in 
the  systems  (not  included  in  report)  some  observations  pertaining  to  the  nature 
of  the  state  of  stress  within  and  the  gross  behavior  of  these  pavements  can  be 
made.  The  transverse  bending  stress  distribution  at  a  distance  from  the  shoulder 
is  linear  with  depth,  with  a  zero  stress  contour  along  the  middle  section. 
Nearer  the  shoulder,  however,  stresses  tend  to  become  uniform  with  thickness, 
primarily  because  the  temperatures  also  tend  to  such  a  uniformity  in  distribution 
from  increased  heat  dissipation  at  the  shoulder  edge.  Because  of  the  anchors  in 
the  anchored  pavement,  stresses  are  locally  increased  and  contours  are  somewhat 
distorted  near  the  anchors  when  compared  with  the  corresponding  conventional 
transverse  bending  stresses. 


126 


fa 
03 
Eh 


a1- 

E-t 


< 

E 

cn 

H 
U 

< 

a: 

H 
6h 

Z 


0.5 


ACTUAL 
IDEALIZED 


A      "I 

s 


rCT 


/ 


/ 


\ 


i 


tt££ 


•A*.*/ 


a  »3    Tsf 
n 


/I 


"^ 


// 


1   in.    =   2.54   cm 

1   TSF   =0.97    kg/ cm' 


f 


a   =1   TSF 
n 


l 
I 

/ 
L 


UNLOADING 


1 


0.05 


0.10 


0.15 


RELATIVE  DEFORMATION ,Ae  (IN.) 
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Longitudinal  stresses  are  quite  different  for  the  two  pavements.  The 
anchors  are  stressed  within  the  soil  and  as  a  result  tend  to  drag  soil  and 
thereby  stress  it  (because  of  frictional  adhesion  at  the  slab/soil  interface). 
The  temperature  difference  between  the  anchors  and  the  soil  caused  by  the  con- 
duction of  heat  between  the  anchors  and  the  slab  make  such  longitudinal  behavior 
possible. 

Thermally  induced  shear  stress  contours  in  the  conventional  pavement 
are  nearly  vertical,  reflecting  the  linear  nature  of  the  transverse  bending 
stresses.  Shear  stress  is  nearly  constant  near  the  shoulder  edge,  especially  at 
night.  The  anchored  pavement  shows  rapid  shear  stress  gradients  near  the  anchors 
due  to  the  restraint  of  curling  deformation  provided  by  the  doweled  connection 
between  the  anchors  and  the  slab. 

To  find  the  effects  of  the  third  dimension  on  the  analysis,  a  single 
slab  was  modeled  (Figures  100  and  101)  with  identical  cross  sections  as  the  cor- 
responding two-dimensional  models.  The  results  for  surface  deflections  of  the 
anchored  pavement  indicate  an  increased  shoulder  (corner  and  edge)  deflection 
(about  double  for  this  particular  condition).  The  interior  of  the  slab  showed  a 
decrease  in  the  deflection  due  to  the  applied  thermal  loads  as  indicated  in 
Figure  102.  The  reduction  in  these  deflections  is  about  50%,  compared  with  the 
two-dimensional  results  in  Figure  98.  Stress  is  expected  to  increase  near  the 
anchors,  the  reduced  deflections  again  being  attributed  to  the  restrictive 
nature  of  these  rigid  anchors. 

The  accuracy  of  the  analysis  would  not  be  credible  unless  the  results 
could  be  checked  against  actual  physical  behavior  of  pavements  in  similar  loading 
environments.  There  are  numerous  road  test  results  published  that  support  the 
results  presented  in  this  study.  In  general,  the  deflections  obtained  agree 
quite  well  for  the  simple  slab,  both  qualitatively  in  general  shape  and  quan- 
titatively in  magnitude.  Differences  can  be  found,  but  it  must  be  remembered 
that  for  a  true  match,  an  identical  physical  or  prototype  system  must  be 
analysed. 

The  results  from  the  analysis  of  the  anchored  pavement  are  obviously 
not  as  easy  to  corroborate  since  there  are  no  prototype  pavements  of  this  con- 
figuration. Because  the  results  obtained  for  the  conventional  pavement  are 
reasonable  and  the  primary  purpose  is  to  study  and  compare  the  two  systems,  the 
results  seem  acceptable. 
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Figure  102.     Deflection  contours  for  three-dimensional 
anchored  slab  during  summer  day.     (in  mils,  down 
is  positive). 
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4.3.6  Moisture 

The  changing  distribution  of  moisture  within  the  load-bearing  layers 
of  pavement  systems  can  often  result  in  variable  behavior  of  the  system  under 
applied  loads.  Mechanical  and  thermal  properties  of  soil  depend  significantly 
upon  moisture  content. 

Some  sources  of  moisture  and  the  impetus  for  movements  are  depicted  in 
Figure  103.  The  more  transient  sources  include  forms  of  seepage  that  result  when 
surface  water  percolates  through  the  surface  layers.  Sources  of  surface  water 
include  precipitation  and  standing  water.  Fluctuations  in  groundwater  table 
elevation  can  change  moisture  contents  in  soil,  and  water  can  reach  upper  layers 
above  a  water  table  by  capillary  action.  Water  vapor  can  also  migrate  through 
partially  saturated  pore  structures  and  reach  load-bearing  layers  in  the 
subgrade. 

There  are  some  trends  in  the  quantity  of  free   moisture  in  any  given 
soil.  The  amount  of  moisture  varies  between  soils  under  similar  conditions. 
Particle  size  and  shape  largely  influence  water  content.  As  particle  size 
decreases,  water  content  generally  increases  due  to  the  larger  specific  surface 
area.  The  thickness  of  film  water  depends  on  the  mineral  structure  and  ions  in 
the  film  fluid.  Particle  surface  charge  can  be  an  important  factor.  Pore 
structure  is  also  important.  Soils  with  fine  pore  structures  are  more  prone  to 
exhibit  capillary  action  than  poorly  graded,  more  uniform  soils.  Other  influenc- 
ing factors  include  the  Atterberg  limits,  viscosity  of  the  permeating  fluid, 
activity,  relative  humidity,  vegetation,  infiltration  quantity,  evaporation,  and 
drainage  conditions. 

Angen  (1973)  has  performed  a  regression  analysis  to  determine  factors 
influencing  moisture  content.  For  road  surface  materials,  moisture  content 
increases  when  density,  the  distance  to  the  groundwater  table,  and  the  percentage 
of  fines  decrease.  For  soils,  moisture  content  increases  as  density  decreases 
and  as  precipitation,  percentage  of  fines,  and  liquidity  index  increase. 

The  moisture  content  of  a  natural  soil  is  altered  when  a  pavement  is 
constructed.  Figure  104  shows  the  variation  in  water  content  of  a  soil  subgrade 
before  and  after  80  days  of  traffic.  The  pounding  of  the  pavement  by  the  dynamic 
vehicle  loads  as  well  as  the  alteration  in  thermal  behavior  of  the  system  create 
the  increase  in  moisture  content  near  the  surface,  in  this  case,  within  3  ft. 
(0.9  m)  of  the  surface. 
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Figure  103.     Sources  of  water  and  directions  of  flow  in 
pavement  systems   (Low  and  Lovell,  1959). 
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Figure  104.     Effect  of  traffic  on  moisture  content  of  subgrade 
soils   (Lytton  and  Kher,  1970). 
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The  behavior  of  soils  and  the  pavement  system  in  particular,  depend 
on  the  water  content  of  subgrade  soil.  Generally  soils  become  weaker  as  water 
content  increases,  especially  in  more  sensitive  clays  where  small  changes  in 
water  content  can  greatly  alter  strength  characteristics.  Because  moisture 
content  in  freeze/thaw  areas  is  greatest  in  the  thaw  period,  the  critical  condi- 
tion to  study  is  the  period  of  rapid  strength  loss  in  spring  with  the  higher 
moisture  content  and  lower  strength.  Modulus  of  the  soil  decreases  under  such 
conditions,  and  Poisson's  ratio  approaches  a  value  of  0.5. 

4.4  ANALYSIS  OF  MECHANICAL  LOADS 

4.4.1  General 

The  primary  effort  in  analytical  investigations  of  pavement  behavior 
has  been  in  the  determination  of  the  effects  of  vehicle  (mechanical)  loads. 
The  deflections  of  the  pavement  surface  are  relatively  easy  to  measure  in  the 
field.  As  techniques  for  measuring  surface  strains  were  developed,  it  became 
possible  to  determine  the  state  of  stress  within  the  pavement,  and  verification 
of  rapidly  developing  theories  for  stress  analysis  of  slabs  was  possible. 

4.4.2  Description  of  Models 

All  models  used  an  identical  18  kip  load  (8100  kg)  in  a  two-wheeled 
axle  configuration.  Each  tire  was  simulated  by  a  single  concentrated  load. 
It  is  realized  that  prototype  loads  are  distributed  over  several  square  inches 
of  area,  but  this  poses  a  problem  in  three-dimensional  analysis.  A  huge  number 
of  additional  elements  would  be  required,  and  the  analysis  would  be  extremely 
costly.  The  single  concentrated  load  does  represent  a  worst  case,  and  stresses 
in  the  immediate  vicinity  of  the  load  will  tend  to  be  higher  than  those  from 
a  distributed  pressure  load. 

Each  three-dimensional  analysis  uses  the  same  layout  of  elements  in 
a  transverse  section  (Figure  105).  An  isometric  view  of  a  100  ft.  (30.5  m) 
slab  (continuous)  model  is  shown  in  Figure  106.  The  bottom  surface  is  fixed 
against  movement  while  vertical  sides  are  on  rollers  such  that  no  out-of -plane 
motion  is  allowed.  These  are  standard  boundary  conditions  used  on  all  succeeding 
three-dimensional  models.  Both  anchored  and  conventional  pavements  are  analyzed 
in  this  configuration.  Additionally,  the  anchored  pavement  has  also  been  compared 
with  conventional  pavement  having  a)  the  same  total  cross  sectional  area  as 
the  anchored  pavement,  and  b)  the  same  rigidity  as  the  anchored  pavement 
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1   ft.    =   30.48   cm 


Figure  106.     Three-dimensional   continous  pavement  model , 
anchored  and  conventional. 
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(Figure  107(a)  through  (e)).  The  effects  of- moisture  variation  are  indirectly 
taken  into  account  by  the  generally  weakening  of  the  top  4  ft.  (122  m)  of  soil  as 
indicated  earlier  and  as  shown  in  Figure  107(e).  As  a  result,  a  two-layer 
elastic  system  is  set  up  in  the  identical  100  ft.  (30.5  m)  slab  configuration 
with  the  same  load.  The  soil  in  this  weakened  zone  is  reduced  to  a  value  of 
100  psi  (7  kg/cm).  With  this  radical  alteration,  the  effects  of  subgrade 
weakening  on  pavement  surface  deflection  response  could  be  readily  seen.  It 
would  also  be  possible  to  determine  roughly  the  percentage  of  load  transferred  to 
the  soil  by  the  anchors. 

Mechanically  uncoupled  jointed  pavements  are  also  considered.  In  this 
analysis  the  load  positions  are:  an  interior  load  and  a  joint  load.  In  addi- 
tion, slab  size  is  considered  with  a  single  12.5  ft.  (3.8  m)  slab  and  a  50  ft. 
(15.3  m)  slab  modeled.  The  single  slab  received  interior  and  joint  loads;  the 
long  slab  (nearly  continuous)  received  only  joint  loads,  which  are  deemed  more 
critical.  The  case  of  interior  load  has  already  been  considered  in  the  100  ft. 
(30.5  m)  slab  analysis. 

Sliding  interface  elements  were  again  employed  for  modeling  slab  con- 
tact with  soil  and  were  used  at  every  node  on  the  subsurface  of  the  slabs.  The 
only  possible  joint  load  transfer  was  through  movement  of  the  subgrade  soil.  A 
jointed  pavement  with  no  load  transfer  across  a  joint  that  is  loaded  at  the  joint 
will  exhibit  somewhat  different  behavior  near  the  joint,  as  depicted  in 
Figure  108.  The  end  surfaces  will  be  unstressed  and  as  a  result  will  rotate  as 
indicated.  Continuous  pavements  have  a  continuity  of  the  middle  plane  beneath 
the  load,  and  since  the  rigidity  is  continuous  in  the  longitudinal  direction, 
will  be  stressed  in  a  plane  passing  through  the  load. 

4.4.3  Results 

Surface  deflections  and  stress  distributions  are  presented  for  various 
conditions  considered.  Figures  109  and  110  present  a  comparison  of  transverse 
and  longitudinal  deflections  of  the  anchored  pavement  with  a  conventional  slab 
(10  in.  (25.4  cm)  thick,  without  anchors),  a  conventinal  slab  with  equivalent 
area  (1.33  x  32  ft.  (0.4  x  9.75  m))  and  a  conventional  slab  with  equivalent 
stiffness  (3  x  32  ft.  (0.9  x  9.75  m)).  The  figures  indicate  the  following 
magnitude  of  deflections  under  load: 


141 


13 


»•  »••• 


a)  Anchored  pavement  (for  dimensions  see  figure  1) 


1  in.  s  2.54  cm 


b)  Conventional  pavement  10  in.   thick. 
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c)  Conventional  pavement  16  in.  thick 
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d)  Conventional  pavement  36  in.  thick 
(equivalent  rigidity  to  a). 
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Figure  107.  Pavement  configurations 
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a)  Continuous 


b)  Jointed. 


Figure  108.  Difference  between  deflection  behavior  with 
load  at  a  joint  in  a  jointed  pavement  and  with  load  on 
a  continuous  pavement. 
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A    CONVENTIONAL  (16 -in.  slab) 

O    CONVENTIONAL  (36 -in.  slab) 

1  mil  =   0.001   in. 
1   in.    ■   2.54   cm 


Figure  109. 


Transverse  surface  deflections  at  load  in 
100-ft.  (30.5-m)  models. 
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a)  anchored  pavement  t2  mils  (1  mil  *  2.54  x  10*3  cm) 

b)  conventional  pavement  29  mils 
(10  In.  *  25.4  cm  slab) 

c)  conventional  with  equivalent      20  mils 
area  (16  in.  »  40.64  cm  slab) 

d)  conventional  with  equivalent        9  mils 
stiffness  (36  in.  *  91.4  cm 

slab) 

The  cross  sectional  concrete  area  for  the  anchored  pavement  is  42.67  ft. 

2  ? 

(3.96  m  )  while  that  of  the  equivalent  stiffness  pavement  is  96  ft.  (about 

9  m  ),  which  1s  more  than  double.  Secondly,  the  36-in.  (91.4  cm)  slab  even  if  it 
is  embedded  24  1n,  (61  cm)  in  subgrade  will  be  heavily  influenced  by  the  freezing 
and  thawing  action  of  the  soil.  The  investigators  strongly  believe  that  advan- 
tages offered  by  the  anchored  pavement  in  terms  of  transferring  the  load  to  below 
freeze/thaw  depth  by  far  excel  the  other  configurations.  In  the  following 
discussions,  therefore,  comparison  has  been  made  only  of  the  anchored  pavement 
with  the  conventional  one  without  the  anchors  (I.e.,  10  1n.  a  25.4  cm  thick). 

Deflection  results  agree  quite  favorably  with  prototypical  behavior. 

Recalling  the  results  presented  by  Moore  et  ah  (1969)  for  a  100-lb.  (484  kg) 

load  and  scaling  for  a  9000-lb.  (4080  kg)  wheel  load,  the  deflection  results  are 

quite  similar  in  magnitude.  (Differences  arise  because  the  soil  and  structural 
properties  are  not  the  same.) 

Surface  deflection  contours  for  the  100-ft.  (30.5  m)  slab  models  are 
shown  1n  Figure  111.  The  conventional  pavement  exhibits  the  typical  concentric 
contour  deflection  basin  often  reported  1n  the  literature.  The  anchored  pave- 
ment, with  its  much  higher  longitudinal  rigidity,  tends  to  "barrel"  and  reveals 
more  cycl1ndr1cal  rather  than  spherical  deflection  patterns.  Deflections  of  the 
anchored  pavement  are  approximately 'one-half  of  the  values  for  the1 conventional 
pavement  at  the  load  points,  but  the  anchored  pavement  deflections  are  somewhat 
larger  near  the  end.  Longitudinal  deflection  of  the  outer  wheel  path 
(Figure  112)  reveals  the  shapes  of  the  deflection  basins  in  the  100-ft.  (30.5  m) 
pavements.  Since  deflections  at  50  ft.  (15.2  m)  from  the  load  were  quite  small 
for  the  conventional  pavement,  it  is  felt  that  boundary  conditions  are  satis- 
factory; deflections  at  the  loads  would  still  be  slightly  less  for  a  truly 
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SYMMETRIC 


LOAD  LOCATIONS 


0.001    in. 
2.54   cm 


Figure  111.     Surface  deflection  contours   (.in  MILS)  for 
continuous  anchored  pavement  [anchor 
locations  are  shown  in  dashed  lines). 
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infinite  continuous  slab.  The  anchored  pavement,  however,  possesses  a  nearly 
uniform  deflection  basin.  Also,  while  the  maximum  surface  deflection  of  the 
anchored  pavement  with  subgrade  weakening  increased  by  15%,  the  maximum  deflec- 
tion of  the  conventional  slab  increased  by  48%.  Thus,  in  the  transference  of 
load  to  lower  layers,  anchors  are  extremely  important,  especially  if  they  are 
found  below  the  frost  line  and  the  zone  of  significant  spring  thaw  weakening. 
Transverse  deflections  of  the  continuous  anchored  pavement  indicate  much  smaller 
deflections  that  are  more  uniform  along  the  transverse  section  (Figures  113 
and  114).  Farther  from  the  load,  conventional  slab  deflections  are  smaller,  but 
shoulder  edges  tend  to  lose  contact  with  the  soil. 

As  expected,  the  deflections  of  the  jointed  pavements  are  substan- 
tially greater  (as  much  as  50%)  than  the  continuous  pavement.  Deflections  are 
also  larger  at  the  joint  for  joint  loads.  The  jointed  pavement  anchor  lengths 
are  the  same  as  those  for  the  continuous  model  and  result  in  similar  total 
stiffnesses.  The  anchors  tend  to  spread  load  over  a  larger  area  and  thus  do  not 
deflect  as  much  as  the  conventional  pavement.  With  an  infinite  (larger)  length 
of  anchors  and  subgrade  extent,  the  edge- loaded  anchored  slab  deflections  should 
be  somewhat  smaller.  The  same  holds  true  for  the  conventional  slab  with  a  larger 
subgrade  extent:  joint  load  deflections  should  be  smaller  at  the  joints  with 
longer  slabs.  In  Figure  115,  the  longitudinal  deflections  of  the  single  slab  are 
shown  with  both  joint  load  and  interior  load.  For  comparison  purposes,  the 
100  ft.  (30.5  m)  results  are  also  included.  Joint  loads  increase  surface  deflec- 
tions in  the  conventional  pavement  significantly  while  they  do  not  affect  the 
anchored  pavement  nearly  so  much.  The  anchors  in  this  case  are  continuous  for 
100  ft.  (30.5  m),  and  the  deflections  are  only  increased  in  the  wheel  path,  not 
directly  over  an  anchor.  Anchor  deflections  are  similar  for  continuous  or 
jointed  slabs.  Had  the  anchors  also  been  discontinuous,  much  larger  deflections 
would  have  occurred,  but  would  be  still  less  than  the  conventional  slab 
deflections. 

The  longitudinal  deflections  of  the  longer  (50  ft.  (15.3  m))  pavements 
are  shown  in  Figure  116.  This  analysis  was  performed  with  the  weakened  subgrade 
condition,  and  very  large  deflections  at  the  joint  in  the  conventional  system 
were  seen.  The  difference  in  the  deflection  of  the  loaded  and  unloaded  slab  is 
apparent,  with  load  transfer  taking  place  via  subgrade  movement.  This  amounts  to 
about  a  50%  transference  of  deflection  but  not  necessarily  of  load.  There  is,  in 
essence,  no  direct  load  transfer.  The  anchored  slab  shows  greater  deflection, 
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(SOLID  SYMBOL  INDICATES 
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Figure  113. 


Transverse  deflections  at  load  in  100-ft. (30.5-ro)  models 
(weak  upgrade). 
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o  ANCHORED 

a  CONVENTIONAL 

(SOLID  SYMBOL  INDICATES 
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1  mil  a  0-.001   in. 
1   in.  s  2.54  an 


Figure  114.     Transverse  deflections  75  in.  from  load 
in  100-ft.  (30.5-m)  models. 
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but  the  anchor  deflection  is  again  identical-  to  that  of  the  continuous  100-ft. 
(30.5  m)  slab  in  the  weakened  soil. 

Tranverse  surface  deflections  are  shown  in  Figures  117  and  118  for 
conventional  slab.  Extreme  deflections  can  be  seen.  The  anchored  pavement 
significantly  decreases  deflections  and  alters  the  nature  of  the  slab  deflec- 
tions so  that  more  uniform  and  controlled  patterns  arise. 

Resulting  stress  contours  for  the  various  load  conditions  in  the 
different  models  have  been  studied,  and  the  results  are  presented  in  the  follow- 
ing pages.  Longitudinal  bending  stresses  again  show  different  behavior  in 
anchored  and  conventional  pavements.  Because  surface  deflections  are  lower,  the 
anchored  pavement  bending  stresses  are  higher.  The  anchors  also  tend  to  pull 
soil  at  about  the  2  to  4  ft.  (0.6  to  1.2  m)  level  to  contribute  to  this  increased 
rigidity.  The  longitudinal  surface  stress  contour  patterns  are  also  different  in 
that  the  continuous  anchored  pavement  has  a  stress  reversal  at  about  20  ft. 
(6.5  m)  from  the  load  while  conventional  pavement  longitudinal  bending  stresses 
tend  to  zero  with  distance  from  the  load.  Transverse  surface  bending  stress 
distributions  for  the  two  concepts  differ  greatly.  Conventional  pavement 
stresses  are  more  regular  and  appear  as  nearby  concentric  contour  lines,  because 
rigidity  is  similar  in  all  directions  far  from  the  pavement  edges. 

The  stresses  introduced  into  the  soil  subgrade  give  further  Insight  on 
the  differences  between  the  two  pavement  concepts.  Vertical  stresses  for  the 
conventional  pavement  appear  as  concentric  rings  at  both  the  24  and  48  1n.  (60 
and  120  cm)  levels.  The  anchored  pavement,  however,  tends  to  distribute  loads 
more  uniformly.  Localized  vertical  stress  Intensification  at  the  24  1n.  (60  cm) 
level  results  from  increased  stress  1n  the  anchors  and  the  assumed  compatibility 
between  anchor  and  soil  that  also  absorbs  some  vertical  stress.  In  the  con- 
tinuous anchored  pavement,  the  lower  level  stress  directly  beneath  the  anchors 
reveals  a  rather  uniform  longitudinal  distribution,  Indicated  by  the  contour 
lines  that  parallel  the  pavement  length.  The  contours  that  cross  the  anchor 
lines  (being  spaced  farther  apart)  indicate  the  more  even  distribution  of  load  1n 
the  soil.  Maximum  vertical  stress  at  24  1n.  and  48  in.  (60  cm  and  120  cm)  depth 
is  also  lower  in  the  anchored  system. 

Jointed  pavement  behavior  reflects  some  similar  trends  in  the  state  of 
stress.  The  jointed  anchored  pavement  tends  to  distribute  loads  deeper  within 
the  subgrade  in  a  more  uniform  manner,  more  uniformly,  in  fact,  than  the  con- 
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Figure  117.  Transverse  deflections  at  joint  in  50- ft.  (15-m) 

pavements  in  weak  soil. 
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Figure  118.  Transverse  deflections  75  in.  from  load  (joint) 
in  50-ft.  (15  m)  jointed  pavements  in  weak  soil. 
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tinuous  anchored  pavement,  as  little  load  is  transmitted  to  adjoining  slabs;  most 
load  is  redistributed  by  the  anchors.  Stresses  are  generally  increased  in  the 
anchored  structural  elements,  but  subgrade  soil  stresses  are  reduced.  Because 
the  load  on  the  conventional  slab  is  applied  as  traveling  in  the  shorter  direc- 
tion, the  longitudinal  bending  stress  distribution  is  no  longer  a  series  of 
simple  concentric  rings.  There  are  also  more  stress  reversals  in  the  jointed 
conventional  slab  than  in  the  continuous  conventional  pavement. 

Vertical  stresses  within  the  subgrade  soils  show  similar  trends.  At 
the  24-in.  (60  cm)  level,  conventional  pavement  stress  contours  are  regular 
rings  while  the  anchored  slab  pattern  indicates  stress  pockets  centered  in  the 
anchors  that  tend  to  pull  the  soil,  spreading  stress  to  the  soil  in  this  local- 
ized area.  Below  the  anchors,  however,  conventional  pavement  stresses  remain 
annular  while  anchored  pavement  stresses  are  nearly  uniform  in  the  longitudinal 
direction,  regardless  of  where  the  load  was  applied  (at  a  joint  or  within  a 
slab). 

4.5  SUPERPOSITION  OF  ENVIRONMENTAL  AND  MECHANICAL  LOADS 

4.5.1  General 

Analysis  of  the  general  response  of  pavements  requires  the  super- 
position of  effects  caused  by  environmental  loads  on  mechanical  (vehicle)  loads. 
A  careful  examination  of  the  predicted  magnitudes  of  the  effects  under  each  type 
of  loading  can  provide  a  better  understanding  of  the  systems  overall  response  and 
subsequent  adequacy  for  designed  use.  Critical  conditions  and  primary  loads  can 
be  ascertained,  and  such  contingencies  can  be  managed.  This  section  seeks  to 
combine  the  effects  of  the  previous  analysis  for  a  more  complete  description  of 
each  system's  gross  response  to  the  types  of  loads  and  loading  combinations  that 
may  be  encountered  throughout  the  design  life  of  the  structure. 

4.5.2  Results 

Of  primary  interest  and  by  far  the  easiest  quantities  to  superimpose 
are  pavement  surface  deflections.  Pavement  analyses  have  historically  begun 
with  prediction  and  examination  of  deflections.  The  results  of  the  superposition 
of  some  extreme  conditions  are  shown  in  Figures  119  and  120.  The  results  clearly 
show  the  significance  of  environmentally  (thermally)  induced  deflections  -- 
curling.  Deflections  from  curling  can  be  as  much  as  one  order  of  magnitude 
larger  than  vehicle- imposed  deflections.  These  strains  are  usually  applied  over 
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Figure  119.  Results  of  surface  deflection  load  superposition  (summer) 
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longer  periods  (several  hours),  whereas  vehicle  loads  are  dynamic  and  nearly 
instantaneous.  Dynamic  loads  are  often  smaller  in  magnitude  than  static  loads 
but  can  cause  reverses  in  deflection  curves.  These  dynamic  vehicle  loads  are 
dissipated  quickly.  Thermally  induced  strains,  however,  do  not  dissipate  so 
quickly  and  can  induce  various  creep  and  fatigue  phenomena,  which,  coupled  with 
constant  vehicle  load  pounding,  can  accelerate  the  deterioration  of  the  system. 

Generally,  the  anchored  pavement  reduces  deflections  by  as  much  as  one 
order  of  magnitude.  Especially  affected  are  curling  deflections  which  are  nearly 
eliminated  by  the  stiffening  of  the  pavement  edges.  When  deflection  is  reduced, 
however,  stresses  are  increased  near  the  anchors,  because  of  the  reaction 
restraint  they  impose  on  the  slab.  It  may  be  necessary  to  include  reinforcing 
rods  as  well  as  steel  fibers  to  adequately  provide  for  such  stress  increases. 

4.6  ANALYSIS  OF  JOINTS 

4.6.1  General 


In  areas  prone  to  large  diurnal  or  seasonal  variations  in  temperature, 

pavement  expansion  and  contraction  become  a  problem.  Joints  must  be  included  in 

the  pavement  system  to  adequately  manage  any  movement  in  the  system  and  still 
provide  a  relatively  smooth  surface. 

The  anchored  pavement  concept  poses  unique  problems  in  joint  design  and 
construction.  Theoretically,  for  a  truly  continuous  pavement,  separation  of  the 
slab  and  anchors  to  promote  sliding  may  be  necessary  near  expansion/contraction 
joints  because  of  the  difference  in  anchor  and  slab  movements.  However,  an 
anchored  pavement  with  a  sleeper  slab  below  the  joint  may  not  require  such  sepa- 
ration, because  deflections  from  thermal  strains  are  smaller  for  shorter  slabs. 

A  prospective  joint  for  the  anchored  pavement  slab  is  presented  in  the 
following  section.  Analysis  was  undertaken  to  determine  the  behavior  of  the 
joint  under  compression  and  expansion  loadings.  Special  construction  details 
are  also  presented. 

4.6.2  Description  of  Joint  and  Mathematical  Model 

The  joint  concept  chosen  is  based  on  a  system  reported  by  Mogilevich  et 
al .  (1976).  The  joint  consists  of  a  set  of  parallel  vertical  plates  (Figure  121) 
coupled  at  the  tops  and  bottoms  by  an  elastomeric  material.  This  material  was 
assumed  to  have  a  modulus  of  500  psi  (35  !<g/cm  )  and  a  Poisson's  ratio  of   0.47, 
typical  of  many  rubberlike  materials.  The  material  is  rather  weak  to  allow 
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Figure  121.  Joint  concept  proposed  for  use  in  anchored  pavement. 
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considerable  flexibility,  and  because  Poissbn's  ratio  is  taken  near  0.5,  defor- 
mations are  at  constant  volume. 

The  joint  is  assumed  to  rest  on  a  sleeper  slab  cast  integrally  with  the 
anchors.  This  inhibits  water  infiltration  and  should  eliminate  problems,  such  as 
pumping,  that  may  occur.  The  lower  surface  of  the  joint  structure  is  friction- 
less  and  can  slide  freely  on  the  sleeper  slab.  The  elastomeric  material  is 
assumed  to  be  bonded  to  the  steel  plates  such  that  no  gaps  or  peeling  occur. 

The  system  was  given  a  2-in.  (5  cm)  compression  and  extension  from  its 
initial  unstressed  state.  This  loading  condition  imposes  a  uniform  deflection  ir 
the  axial  direction  (no  rotation  or  curling  of  the  slab  end)  that  can  be  con- 
trolled by  the  sleeper  slab  design. 

4.6.3  Results 

Deflections  and  stresses  were  determined  for  a  plane  strain  longi- 
tudinal section  of  a  typical  transverse  joint.  Figure  122  presents  deflections 
resulting  from  the  expansion  and  contraction  loadings.  In  this  figure,  vertical 
deflections  are  exaggerated  by  a  factor  of  2.5  to  1.  The  resulting  deflections 
create  small  ridges  approximately  0.2  in.  (5  mm)  high  and  depressions  of  about 
the  same  magnitude.  This  is  consistent  with  values  reported  by  Mogilevich  (1976] 
for  a  30%  reduction  in  joint  width,  about  0.2  in.  (4.8  mm)  maximum. 

Axial  stresses  in  the  concrete  slab  near  the  joint  induced  by  the 
expansion  and  contraction  of  the  joint  are  also  shown  in  Figure  122.  Maximum 
stresses  of  about  120  psi  (8.4  kg/an  )  in  both  tension  and  compression  are 
manageable  in  steel  fiber  reinforced  concrete  (SFRC).  Care  must  be  taken  to 
ensure  proper  anchoring  of  the  joint  to  the  slab  by  hooks  or  other  conventional 
anchors.  The  pressure  distribution  shape  is  consistent  with  its  loading  geo- 
metry, as  the  low  stress  portion  occurs  in  the  void  area  that  transfers  only  a 
small  portion  of  the  load  through  the  end  steel  plates. 

The  results  of  this  analysis  indicate  the  feasibility  of  this  joint- 
filling  concept.  Some  physical  details  would  need  refining  for  field  testing  anc 
evaluation. 
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CHAPTER  V 
CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  CONCLUSIONS 


This  study  has  attempted  to  further  advance  the  methodology  for  anal- 
ysis and  design  of  pavements.  Of  primary  concern  was  the  effect  of  environmental 
loads  placed  upon  the  system,  in  particular,  cyclic  temperature  differentials. 
Relative  magnitudes  of  the  analytically  obtained  deflections  compared  favorably 
with  field  data  collected  from  various  road  tests  and  reported  in  the  literature. 
It  is  hoped  that  the  use  of  three-dimensional  finite  element  analysis  will  gain 
increasing  acceptance  as  more  efficient  techniques  are  developed  to  handle  the 
associated  large  system  of  equations,  and  also  as  faster  and  more  powerful 
computers  are  utilized  to  cut  analysis  costs.  Many  problems  in  pavement  system 
analysis,  as  well  as  in  geotechnical  engineering  in  general,  can  only  be  solved 
adequately  by  a  full  three-dimensional  analysis. 

The  results  of  the  analysis  indicate  the  significance  of  environ- 
mentally induced  stress.  Pavement  warping  deflection  can  be  greater  than 
deflection  due  to  applied  wheel  loads. 

Various  combinations  of  loading  conditions  that  include  as  components 
the  environmental  and  mechanical  loads  are  possible.  It  is  possible  to  choose  a 
case  that  is  sufficiently  detrimental  and  occurs  often  enough  to  be  considered  a 
worst  case.  Usually  such  a  condition  involves  high  stress  concentrations  and 
numerous  stress  reversals.  A  conventional  slab  that  warps  down  (during  daylight 
hours)  then  is  loaded  mechanically  will  experience  deflections  similar  to  those 
in  Figures  119  and  120,  and  reversals  of  transverse  bending  stresses  will  result 
in  a  "wrinkled"  slab.  The  anchored  pavement  deflections  are  relatively  uniform 
and  smaller,  thus,  will  be  much  more  predictable;  however,  streses  near  the 
anchors  in  the  structural  elements  of  the  anchored  pavement  are  increased. 
Precautions  must  be  taken  to  ensure  adequate  reinforcement  in  these  areas. 

The  anchored  slab  offers  two  distinct  advantages  over  a  conventional 
slab.  First,  deflections  are  lower  and  more  uniform.  Second,  stresses  in  the 
soil  are  lower  and  distributed  more  widely  by  the  rigid  anchors.  A  significant 
portion  of  pavement  distress  mechanism  arises  from  subgrade  failure.  If  loads 
are  transmitted  deeper  within  the  subgrade,  to  soil  that  is  under  greater 
confining  stresses  (and  as  a  result  is  stronger)  and  where  moisture  and 
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temperature  fluctuations  are  not  acute,  subgrade-related  failure  is  less  likely 
to  occur. 

As  with  any  engineering  solution,  there  are  trade-offs.  Generally 
soil  behavior  is  much  more  random  and  uncontrollable;  this  should  be  a  serious 
consideration.  The  actual  structural  elements  can  be  strengthened  in  a  con- 
trollable manner,  and  tolerance  to  increased  stresses  can  be  managed. 

5.2  RECOMMENDATIONS 

Finite  element  modeling  techniques  for  the  analytical  determination  of 
a  pavement  system  response  to  various  loading  conditions  have  been  presented. 
Full  three-dimensional  analysis  of  continuous  and  jointed  pavements  is  rela- 
tively costly,  but  the  analysis  is  much  more  sophisticated  as  complex  loading  and 
boundary  conditions  and  system  geometry  are  more  accurately  modeled.  Plane 
strain  analysis  of  longitudinal  or  transverse  sections  completely  ignores  out- 
of-plane  rigidity  of  the  structure,  which  is  extremely  important,  especially  in 
the  anchored  pavement.  Loads  are  also  input  as  line  loads  of  infinite  extent  and 
do  not  accurately  portray  prototypical  configurations.  Techniques  to  extra- 
polate three-dimensional  behavior  from  a  two-dimensional  model  using  prismatic 
solid  finite  elements  (Cheung,  1976;  Pichumani,  1973)  have  been  developed. 
Continuous  pavements  with  three-dimensional  arbitrary  loading  patterns  approx- 
imated by  Fourier  expansions  in  the  longitudinal  direction  can  be  modeled; 
however,  jointed  or  otherwise  discontinuous  (cracked)  pavements  have  not  been 
adequately  simulated. 

This  study  has  considered  fully  continuous  and  discontinuous  pavement 
slabs  to  simulate  continuous  rigidity  and  load  transfer  and  finite  rigidity  and 
no  load  transfer,  respectively.  These  represent  the  extremes  in  actual  systems; 
the  prototypes  are  actually  categorized  between  continuous  and  discontinuous. 
It  is  possible  to  include  the  effects  of  partial  load  transfer  by  several 
methods.  Perhaps  the  simplest  method  that  does  not  add  any  elements  to  the 
system  involves  simulating  load  (force)  transfer  by  constraining  the  movement  of 
nodes  along  the  joint  of  unloaded  slabs  adjacent  to  a  loaded  slab.  It  is 
possible  to  write  an  equation  in  the  computer  program  utilized  (ANSYS)  to  accom- 
plish this: 


c  a  I  s  a,  Ui  11 


165 


where  c  is  a  constant,  a.  is  factors  and  in"  is  displacements.  For  example,  if  a 
node  2  was  allowed  to  move  in  the  x  direction  only  half  as  much  as  a  node  1,  the 
equation  for  this  constraint  would  be: 

x]>  -  2x2  =  0  12 

so  that  a,  =  1,  a2  -   -2,  and  c  =  0.  It  is  possible  to  simulate  more  complex 
phenomena  such  as  limited  rotational  (moment)  transfer  by  more  complex 
equations. 

Another  method  to  simulate  load  transfer  is  to  include  dowel  elements 
in  the  mesh,  either  in  single  or  double  layers.  Care  must  be  taken  when  using 
truss  elements  with  pinned  end  conditions  to  model  dowels  embedded  in  concrete. 
A  system  of  springs  can  be  used  to  effectively  fix  the  rotation  of  the  bar  with 
respect  to  the  slab  end. 

Other  techniques  can  involve  the  use  of  sufficiently  weak  joint  filler 
material  to  transfer  some  load  and  absorb  the  rest  as  deformations,  networks  of 
springs  to  transfer  load,  or  fractional  loads  on  each  side  of  the  joint. 

The  heat  transfer  analysis  performed  for  this  study  is  a  sufficiently 
accurate  appraisal  of  conditions  above  freezing.  However,  there  are  problems 
associated  with  thermal  properties  of  materials  that  phase  (water  in  soil  that 
freezes  and  thaws),  in  particular  with  the  heat  capacity.  As  stated  earlier,  as 
water  freezes,  latent  heat  is  released,  resulting  in  an  abrupt  discontinuity  in 
the  heat  capacity  -  temperature  relation,  usually  over  a  yery   small  temperature 
range.  Comini  et  al_.  (1974)  present  a  method  for  resolving  this  problem  when 
using  the  finite  element  method  for  nonlinear  heat  conduction  problems.  By 
redefining  specific  heat,  c,  as  the  derivative  of  enthalpy,  H,  such  that: 

c  =  dH  13 

c   dT  13 

and  noting  that  the  enthalpy-temperature  relation  in  the  interval  of  freezing  and 
thawing  temperatures  is  sufficiently  smooth,  heat  capacity  can  be  approximated 
by  determining  the  gradient  of  enthalpy  with  respect  to  temperature.   In  this  way 
any  peaks  in  heat  capacity  relationship  can  be  included. 

It  is  also  possible  to  model  the  complex  phenomenon  of  ice  lense 
formation,  with  specific  attention  focused  upon  the  effects  of  such  formations. 
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Figure  123  presents  a  conceptual  finite  element  mesh  to  model  isolated  lenses  in 
two  or  three  dimensions.  As  water  changes  its  phase  from  liquid  to  solid,  its 
volume  increases  by  about  S%.     As  a  result,  pressures  develop  that  tend  to  push 
pavements  upward.  This  movement  can  be  simulated  by  opening  a  slit  in  the  mesh 
with  double  nodes,  then  moving  the  nodes  on  either  side,  as  shown.  The  amount  of 
movement  depends  on  the  thickness  of  the  ice  Tense,  which,  for  any  particular 
soil,  can  be  determined  experimentally  or  by  field  tests.  For  soils  that  exhibit 
rhythmic  banding  of  ice  in  discrete  layers,  this  is  an  approximate  behavioral 
model.  Figure  124  shows  the  deflection  in  a  conventional  slab  at  a  joint.  For 
the  case  of  an  anchored  slab  an  analysis  was  performed  (in  2  dimensions  only) 
using  the  lense  as  shown  in  Figure  125,  which  also  shows  the  resulting  upward 
movement.  The  pavement  will  probably  exhibit  the  same  behavior  if  the  total 
upward  displacement  of  the  soil  surface  is  the  same. 

Finally,  it  is  also  possible  to  analyze  the  dynamic  response  of 
pavements  due  to  moving  load  by  the  finite  element  method.  Lewis  and  Harr  (1969) 
studied  moving  loads  on  slabs  within  two  dimensions  with  complex  boundary 
conditions,  such  as  slab  warping,  variable  subgrade  modulus,  and  slab  discon- 
tinuities. The  technique  involves  adding  an  acceleration  term  to  the  left  side 
of  the  plate  on  the  elastic  foundation  equation  that  has  the  following  form: 

where  p  and  H  are  density  and  thickness  of  the  slab,  respectively,  and  w  is 
vertical  movement.  A  suitable  finite  difference  formula  is  written,  and  the 
displacement  or  stress  history  of  points  within  the  slab  is  possible. 

The  finite  element  method  can  also  be  used  for  the  transient  dynamic 
analysis  of  the  more  complex  anchored  slab.  Loads  are  applied  at  nodal  points  on 
the  surface  in  ramplike  fashion,  as  shown  in  Figure  126.  With  sufficient  overlap 
of  load  and  spacing  in  time  governed  by  vehicle  velocity,  dynamic  wheel  loads  can 
be  simulated.  A  transient  dynamic  analysis  similar  to  the  transient  heat 
transfer  analysis  is  possible  with  ANSYS. 

It  has  been  reaffirmed  that  the  finite  element  method,  in  either  two  or 
three  dimensions,  can  be  an  extremely  versatile  and  effective  method  of  analysis 
for  many  phases  of  pavement  design,  especially  for  the  critical  environmental 
load  regime.  Its  widespread  use  by  state  highway  design  departments  is 
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Figure  123.  Possible  technique  for  the  analysis  of  the 
effects  of  isolated  ice  lenses  in  two  and 
three  dimensions  by  the  finite  element  method. 
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Figure  125.  Results  of  ice  lense  growth  24  in.  (61  cm)  beneath  pavement 
(anchored  slab). 
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Figure  126.  Scheme  to  study  dynamic  loads  in  pavement 
systems  using  transient  dynamic  finite  element 

analysis. 
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anticipated  as  more  advanced  techniques  are  developed  to  analyze  contingencies 
that  have  traditionally  been  overlooked,  ignored,  or  confined  to  field  or 
laboratory  investigation.  Versatile  yet  usable  codes  can  be  written  with 
sophisticated  data  generation  schemes,  making  model  development  much  quicker  and 
less  costly.  Comprehensive  laboratory  studies  to  determine  realistic  input 
material  properties  and  geometric  seviceability  requirements  coupled  with  engi- 
neering judgement  will  always  be  required  to  design  or  renovate  pavements. 
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APPENDIX  A:   TABLES  OF  SURFACE  STRAINS 


Table  5   Surface  strains  (microin./in.)  of  AS62  far  Load  position  1 


Foil  gage 
number 


Load  (lb.)   (1  lb.  =  0.45  kg) 


250 


500 


750 


100Q 


1250 


1500 


55 

115 

150 

190 

243 

290 

2 

5 

-36 

-50 

-63 

-50 

-50 

3 

30 

57 

40 

46 

50 

CO 

4 

-7 

-14 

-18 

-23 

-25 

-35 

5 

-12 

-22 

-37 

-45 

-54 

-58 

5 

14 

20 

27 

33 

38 

45 

7 

11 

21 

28 

42 

49 

45 

3 

-3 

-4 

5 

11 

13 

12 

9 

4 

-9 

-18 

-19 

-24 

-32 

10 

5 

7 

10 

13 

16 

11 

Table  6  Surface  strains  (microin./in.)  of  AS62  for  load  position  2 

Load  (lb.)   (l  lb.  =  0.45  kg) 


Foil  gage 
number 

250 

1 

10 

2 

10 

3 

-16 

4 

-10 

5 

-10 

6 

-12 

7 

-5 

3 

-9 

9 

-7 

10 

-6 

500 


750 


20 

-32 

IS 

78 

-18 

36 

-13 

-15 

-18  . 

-23 

-17 

-40 

7 

-12 

-16 

-14 

-16 

-21 

-9 

-14 

1000 

-33 
92 

45 
-19 
-25 
-36 
-13 
-24 
-34 
-23 


1250 

-20 
100 
63 
-26 
-30 
-42 
-31 
-20 
-36 
-24 


1500 

15 
126 

66 
-24 
-34 
-33 

-6 
-24 
-40 
-30 
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Table  7   Surface  strains  (microin./in. }  of  AS62  for  load  Dosition  3 

Load  (lb.)    1  lb  =  0.45  kg 

750      1000      1250      1500 

-68  -88       -30       -70 

-84      -100      -HO      -170 


Foil  gage 
number 

250 

500 

1 

-12 

-46 

2 

-48 

-72 

3 

8 

-30 

4 

7 

10 

5 

-4 

-9 

6 

14 

8 

7 

3 

4 

3 

8 

13 

9 

-13 

-24 

10 

-22 

-40 

-37       -46       -25 


o 


13  12  12  10 

-15  -20  -16  -13 

12  15  23  20 

7  9  12  10 

17  20  23  30 

-36  -47  -55  -62 

-64  -83  -100  -116 


"Table  8    Surface  strains  (microin/in. )  of  AS62  for  load  position  4, 


Foil  gage 
number 

Load 

(lb.)      1  lb 

=  0.45  kg 

250 

500 

750 

1000 

1250 

1500 

1 

-230 

-234 

-166 

-149 

-125 

-94 

2 

-227 

-246 

-254 

-269 

-257 

-278 

3 

-260 

-248 

-246 

-259 

-259 

-249 

4 

-53 

-60 

-61 

-66 

-72 

-77 

5 

7 

-38 

-35 

-23 

-12 

3 

6 

-2 

-142 

-138 

-145 

-154 

-150 

7 

-8 

-139 

-145 

-152 

-167 

-182 

8 

-2 

-50 

-54 

-57 

-60 

-60 

9 

-1 

-52 

-51 

-50 

-49 

-44 

10 

-2 

-47 

-46 

-42 

-40 

-39 
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Table  9"    Surface  strains  (microin./in.)  of  AS62  for  load  position  5, 


Foil  gage 
number 

Load  (lb.) 

1  lb  > 

•  0.45  kg 

250 

500 

750 

1000 

1250 

1500 

1 

-181 

-178 

-153 

-127 

-119 

-100 

2 

-204 

-190 

-190 

-168 

-151 

-135 

3 

-276 

-274 

-294 

-288 

-276 

-280 

4 

-52 

-56 

-60 

-64 

-69 

-72 

5 

-48 

-42 

-40 

-35 

-29 

-22 

6 

-144 

-148 

-152 

-154 

-156 

-149 

7 

-155 

-164 

-166 

-169 

-174 

-178 

8 

-46 

-48 

-54 

-56 

-59 

-62 

9 

-42 

-34 

-31 

-27 

-13 

-5 

10 

-40 

-42 

-45 

-48 

-44 

-43 

Table  10    Surface  strains  (microin./in.)  of  AS62  for  load  position  6, 


Foil  gage 
number 

Load  (ir 

i.)   1  lb  * 

0.45  kg 

250 

500 

750 

1000 

1250 

1500 

1 

-23 

-49 

-97 

-112 

-143 

-158 

2 

-35 

-30 

-65 

-63 

-30 

-96 

3 

-28 

40 

35 

32 

38 

42 

4 

3 

6 

12 

10 

16 

14 

5 

3 

4 

3 

2 

-4 

-2 

6 

5 

-17 

-14 

-16 

-14 

-16 

7 

-10 

-13 

-14 

-17 

-16 

-18 

8 

8 

6 

8 

10 

14 

17 

9 

2 

-6 

2 

12 

16 

20 

10 

8 

10 

19 

36 

51 

70 
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e  11    Surface 

strains 

(micro in 

./in. 

)  of 

S62 

fo 

r  load  position  1 . 

Foil  gage 
number 

250 

Loac 

[  (lb 
500 

) 

1  lb  =  0.45  kg 
750 

1 

2 

6 

12 

17 

3 

0 

-1 

-2 

4 

-4 

-5 

-6 

5 

-1 

-2 

-2 

6 

4 

2 

0 

7 

51 

91 

131 

8 

-10 

-19 

-29 

9 

0 

2 

2 

10 

2 

4 

6 

able  12 

Surface 

strains 

(microin 

./in. 

)  of 

S62 

for 

loac 

Foil  gage 
number 

250 

Loac 

1  (lb 
500 

•  ) 

1  lb 

=  0 
750 

1 
2 

-8 

-13 

-17 

3 

-15 

-34 

-53 

4 

-10 

-17 

-21 

5 

-8 

-15 

-21 

6 

8 

10 

10 

7 

-10 

-14 

-10 

8 

6 

10 

12 

9 

15 

34 

51 

10 

6 

10 

12 

182 


Table  13    Surface  strains  (microin./in.)  of  S62  for  load  position  3 


Foil  gage 
number 

Load  (lb.) 

1  lb  * 

250 

500 

750 

1 

2 

-8 

-13 

-14 

3 

16 

26 

35 

4 

-26 

-41 

-49 

5 

-30 

-130 

-150 

6 

2 

8 

8 

7 

-19 

-27 

-23 

8 

8 

12 

14 

9 

-16 

-24 

-32 

10 

21 

34 

39 

Table  14    Surface  strains  (microin./in.)  of  S62  for  load  oosition  4, 

Foil  gage  Load  (lb.)  1  lb  =  0.45  kg 
number 

250  500  750 

1 

2  28  41  58 

3  -6  -8  -10 

4  4  4  4 

5  111 

6  -6  -10  -17 

7  45  78  116 
3  -24  -36  -49 
9  4  8  13 

10  -4  -3  -1 
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Table  15 

Surface 

strains  (micro in 

./in.)  of  S62 

for 

loa 

Foil  gage 
number 

Load  (lb.) 

1  lb 

■  C 

* 

250 

500 

750 

1 
2 

9 

13 

13 

3 

-19 

-42 

-72 

4 

9 

19 

25 

5 

-3 

-5 

-11 

6 

2 

1 

'  0 

7 

5 

15 

23 

8 

-9 

-13 

-19 

9 

17 

40 

64 

10 

-11 

-19 

-30 

Table  16    Surface 

strains  (microin 

./in. 

)  of  S62 

for  loa 

Foil  gage 
number 

Load  (lb.) 

1  lb  = 

250 

500 

750 

1 

2 

-3 

-9 

-13 

3 

16 

26 

36 

4 

-11 

-16 

-21 

5 

-22 

-28 

-17 

6 

-1 

1 

5 

7 

-13 

-19 

-26 

8 

3 

6 

11 

9 

-15 

-25 

-34 

10 

-7 

11 

14 
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Table  17    Surface  strains  (microin./in.)  of  ASE52  for  load  position  1 

Foil  gage  Load  (lb.)    1  lb  -  0.45  kg 
number 

100         200  300 

1  23  54           90 

2  -2-2  2 

3  1  3           12 

4  -4  -3  -15 
5-123 

6  -1  -1           -2 

7  -2  -2  -4 
3  -2  -2  -5 
9  Oil 

10  0           0-1 


.Table  18    Surface  strains  (microin./in.)  of  ASE62  for  load  position  2. 

Foil  gage 
number 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 


Load   (lb.) 

1  lb  =  0.45  kg 

00 

200 

300 

3 

-2 

5 

2 

8 

17 

-1 

-6 

-2 

0 

0 

-2 

0 

0 

0 

2 

0 

1 

2 

2 

2 

0 

-2 

-4 

0 

-2 

0 

0 

0 

2 
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Table  19    Surface 

strains  (micro in 

./in. 

)  of  ASE52 

for 

loa 

Foil  gage 
number 

Load  (lb.) 

1 

lb 

»  0. 

1Q0 

200 

300 

1 

*10 

-21 

-54 

2 

-a 

-19 

-30 

3 

-2 

-4 

-14 

4 

2 

6 

12 

5  • 

-2 

0 

1 

6 

0 

2 

4 

7 

-1 

0 

0 

8 

2 

6 

10 

9 

-2 

-4 

-4 

10 

-4 

-6 

-8 

Table  20    Surface  strains  (microin./in. )  of  ASE52  for  load  position  4. 


Foil  gage 
number 

Load  (lb.)   1 

lb  =  0.45 

10O 

200 

300 

1 

-15 

-21 

-13 

2 

0 

3 

6 

3 

2 

8 

15 

4 

-3 

-6 

-10 

5 

4 

12 

19 

6 

-1 

-2 

-4 

7 

-4 

-10 

-14 

8 

-2 

-3 

-6 

9 

2 

6 

6 

10 

0 

2 

0 
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Table  21    Surface  strains  (microin./in. )  of  ASE62  for  load  position  5 


Foil  gage 
number 

Load  (lb.) 

1  lb  = 

100 

200 

300 

1 

-a 

-10 

-15 

2 

4 

8 

15 

3 

-2 

-4 

0 

4 

-2 

-4 

-7 

5 

2 

5 

9 

6 

0 

-2 

-4 

7 

-2 

-6 

-9 

8 

-2 

-4 

-6 

9 

4 

6 

11 

10 

0 

2 

2 

Table  22    Surface 

strains  (micro in 

./in.)  of  AS 

£52 

for  loa 

Foil  gage 
number 

Load  (lb.) 

1 

lb  =  0. 

100 

200 

300 

1 

-23 

-41 

-77 

2 

-4 

-tlO 

-20 

3 

0 

6 

1 

4 

0 

2 

4 

5 

0 

2 

2 

6 

2 

2 

4 

7 

0 

-2 

-4 

8 

2 

4 

8 

9 

4 

9 

14 

10 

10 

21 

30 
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APPEMOIX  8:  GEOTECHNICAl  PROPERTIES  OF  SU3GRADE 

The  properties  of  the  subgrade  used  in  experimental  investigations 
were  determined  by  a  plate-bearing  test,  unconfined  compression  tests  (UC), 
triaxial  compression  tests,  consolidation  tests,  and  soil  classification  and 
subgrade  behavior. 

PLATE-8EARIN6  TEST 

The  soil  subgrade  was  loaded  by  a  round  rigid  plate.  Observed  settle- 
ment and  consequent  strains  were  plotted  against  theoretically  computed  stresses 
to  estimate  the  modulus  of  elasticity  of  soil. 

The  sampling  area,  as  shown  in  Figure  4  of  the  steel  container,  was 
used  to  perform  the  Plate  Bearing  Test. 

The  dimensions  of  the  sampling  area  are  36.25  in.  (92  an)  wide, 
23.38  in.  (61  an)  long,  and  23.75  in.  (60  cm)  deep.  Eleven  layers  of  silty  clay 
soil  were  placed  and  compacted  in  the  same  way  and  simultaneously  with  the 
testing  area,  so  that  similar  properties  can  be  expected  of  this  subgrade.  The 
height  of  soil  subgrade  was  measured  to  be  22.75  in.  (57.8  cm). 

A  pair  of  sensors  (1-in.  (2.54  cm)  diameter)  were  placed  in  the  soil  at 
a  depth  of  1.63  in.  (4.14  cm)  in  the  center  of  the  sampling  area.  The  average 
distance  between  the  two  sensors  was  found  to  be  about  0.75  in.  (1.9  cm).  Prior 
to  the  placement  of  the  sensors,  a  standard  calibration  procedure  was  followed. 
A  soil  strain  gage  (Figures  127  and  128)  and  a  calibration  fixture  for  sensors 
were  used. 

Method  of  Calibration 

One  sensor  was  placed  in  the  sliding  rod  and  the  other  in  the  coarse 
adjustment  of  the  calibration  fixture  as  showr  in  Figure  129.   They  were  con- 
nected in  the  signal  input  and  in  the  oscillation  output  of  the  soil  strain  gage 
respectively.  The  calibration  signal  was  set  to  14,  the  coil  separation  to  1 
diameter,  and  the  sensitivity  at  low  level  (2  turns  clockwise).  The   calibration 
was  performed  without  balancing  the  phase. 

The  initial  distance  between  the  two  sensors  was  measured  and  the 
corresponding  amplitude  was  recorded.  A  change  in  the  micrometer  corresponded  to 
an  amplitude  change  which  set  the  dial  meter  to  zero.  Figure  130  shows  the 
relationship  between  the  sensor  spacing  and  the  amplitude  reading. 

After  the  top  soil  surface  was  leveled,  a  rigid  steel  footing  was 
placed  in  the  center  of  the  sampling,  area.  The  dimensions  of  the  footing  were 
5.5-in.  (14  cm)  diameter  and  2  in.  (5  on)  thick.  Three  dial  gages  were  attached 
from  the  sides  of  the  container  at  120  degree  angles,  and , they  were  touching  the 
footing  to  measure  the  surface  deflections. 

An  average  deflection  was  computed  from  the  three  measured  deflec- 
tions. A  point  load  was  applied  at  the  center  of  the  footing.  The  point  load  was 
connected  to  a  manual  jack  through  a  load  cell.  The  load  cell  was  connected  to  a 
power  supply  and  to  a  multimeter.  The  net  point  load  applied  on  the  footing  is 
computed  as  follows: 
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Net  load  (lb.)  =  Multimeter  reading  (uV.)  +  lg  ^ 

To  find  the  stresses  in  the  soil  at  1.63-in.  (4.14  cm)  depth  underneath 
the  load,  the  theory  of  circular  load  on  semi-infinite  mass  was  used.  The 
following  notations  were  considered: 

o       3    vertical  stress  at  1.63  in.  (4.14  cm)  depth 

z       3    distance  between  top  surface  and  point  where  vertical 
stress  is  to  be  computed  *  1.63  in.  (4.14  cm) 

a  3  radius  of  footing  -   2.75  in.  (7  cm) 

Pav  s  average  applied  pressure 

P  -  total  load 

e  s  eccentricity  and  is  zero  in  our  case 

p  s  vertical  surface  displacement 

£  s  secant  modulus  of  soil 

v  -  Poisson's  ratio  of  soil 

C       3    |3  0.59. 
a 

From  Poulos  and  Davis   (1974)   and  noting  that  e  s  0  and  5  =  0.59, 
we  conclude  that: 

£2L  *  0.56  B"2 

pav 

Note  that: 

.  _p  B-3 

*a< 


Hav  2 


From  equations   (B-2)   and  (3-3): 

az  (psi)   -  0.02357  P  (lb..).  B-4 

Equation  (B-4)  is  a  relationship  between  a  point  load  P  and  the  cor- 
responding stress  a    at  1.63-in.  (4.14  cm)  depth.  The  percentage  strain  at  this 
depth  was  defined  as  the  difference  of  the  original  sensor  spacing  with  no  load 
from  the  one  at  a  given  load,  divided  by  the  original  one,  and  the  result 
multiplied  by  100. 

Figure  131  shows  a  plot  of  applied  load  versus  average  surface  deflec- 
tion for  ten  load  cycles. 
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Figure  132  shows  a  plot  of  vertical  stress  versus  average  strain  at 
1.63  in.  (4.14  an)  depth  for  10  load  cycles.  After  the  tenth  load  cycle,  the 
footing  was  loaded  up  to  failure  (Figure  133) .  The  stress  at  failure  was  20  psi 
(1.4  kg/cm  )  and  the  corresponding  strain'was  18%.  After  unloading,  the  residual 
plastic  deformation  was  observed  as  0.1  in. /in.  (an/cm). 

Using  the  slope  of  the  tenth  cycle  of  the  stress-strain  curve-  the 
secant  modulus  of  the  subgrade  comes  out  to  be  around  850  psi  (60  kg/ cm  ) 
(Figure  134.) .  The  relationship  between  the  vertical  surface  deflection  and  the 
secant  modulus  of  the  soil  is  quoted  as: 

fl  .Tm   2,   pav  a  8-5 

Equation  (8-5)  becomes: 

*z  .  f  (1  -  v2)  ^|_1  8-6 

Equation  3-3  becomes: 

Ap  B-7 


^v'^I 


ira 
Ap     and  Ap  were  found  from  the  tenth  load  cycle  in  Figure  131    as: 

Ap2      =      0.066  in.   (0.17  cm) 

Ap  *      284  lbs.   (128  kg) 

Therefore,  by  combining  equations  8-6  and  B-7,  the  result  is: 

E  (psi)   =  782.37   (1  -  v2)  B-8 

Secant  modulus  is  plotted  against  Poisson's  ratio  in  Figure  135.  It 
may  be  noted  that  for  v  s  0.45,  the  secant  modulus  is  623  psi  (44  kg/cm*). 

UNCONFINED  COMPRESSION  TEST  (UC) 

Five  unconfined  compression  tests  (ASTW-02166)  were  completed  to  eval 
uate  the  compressive  strength  of  soil.  Many  empirical  relationships  exist  to 
relate  E  with  undrained  unconfined  compressive  strength. 

The  unconfined  compression  test  is  the  simplest  and  quickest  way  to 
measure  the  unconfined  compressive  strength  of  a  cylinder  of  soil.  No  lateral 
support  is  offered.;  therefore,  it  is  used  only  for  cohesive  soils; 

The  failure  surface  tends  to  develop  in  the  weakest  portion  of  the 
sample.  Internal  soil  conditions  such  as  degree  of  saturation  and  pore  water 
pressure,  cannot  be  controlled.  The  sample  is  unprotected  from  moisture  lost 
during  testing.  The  sample  was  extruded  from  the  sampling  area  by  a  modified 
Harvard  compaction  unit.  The  procedure  followed  the  ASTM  specifications.  The 
samples  had  an  initial  diameter  of  1.36  in.  (3.45  cm)  and  an  initial  height  of 
2.75  in.  (6.96  cm). 
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A  plot  of  the  applied  stress  versus  the  percentage  average  strain  is 
given  in  Figure  136.  The  secant  modulus  versus  the  percentage  average  strain 
was  also  plotted  (Figure  137). 

-The  maximum  value  of  the  secant  modulus  is  observed  to  be  25  psi 
(18  kg/an)  at  a  strain  of  0.0184  in. /in.  (on/cm).  At  0.1%  strain,  the  cor- 
responding secant  modulus  is  about  155  psi  (11  kg/cm  ). 

TRIAXIAL  TESTS  (CU) 

The  tests  that  were  done  are  called  Q  or  consol idatad-undrained  or 
consol i dated-quick. 

The  procedure  followed  was  as  per  ASTM  specifications. 

The  samples  used  were  about  1.3  in.  (3.3  cm)  in  diameter  and  2.8  in. 
(7.1  cm)  in  height.  The  chamber  pressures  applied  were  1,  5,  25,  50,  75,  and 
100  psi  (1  psi  s  0.07  kg/an  ).  Two  sets  of  tri axial  tests  were  performed, 
depending  on  the  way  the  samples  were  taken.  Both  tests  measured  the  compressive 
strength  of  a  cylinder  of  soil  at  different  lateral  pressures. 

CU  Test  on  Samples  Extended  From  the  Sampling  Area  by  a 
Modified  Harvard  Compaction  Mould  (Non-Standard)" 

_.  A  plot  of  the  deviator  stress  versus  the  percentage  strain  is  given  in 
Figure  138.   Figure  139  shows  the  Mohr  circles  at  failure.  The  angle  of 
friction  is  12°  and  the  cohesion  is  500  psf  (0.24  kg/cm  ).  The  modulus  of 
elasticity  of  the  soil  versus  the  chamber  pressure  is  shown  in  Figure  140.   The 
modulus  of  elasticity  of  0  psi  confining  pressure  is  taken  as  the  secant  modulus 
from  the  UC  test  at  0.1%  strain.  The  rest  points  were  taken  on  the  initial 
tangent  values  of  stress-strain  curve.  At  about  2  psi  (0.14  kg/cm  )  confining 
pressure,  the  corresponding  modulus  of  elasticity  of  the  soil  is  450  psi 
(31.5  kg/a/). 

CU  Test  on  Undisturbed  Samples  Taken  From  the  Testing  Area 

A  plot  of  the  deviator  stress  versus  percentage  strain  is  given  in 
Figure  141.,  Figure  142  shows  the  Mohr  circles  at  failure.  The  angle  of 
friction  is  18°  and  the  cohesion  is  200  psf  (0.10  kg/cm  ).  The  secant  modulus  of 
the  soil  was  defined  as: 

E  .  *d  B-9 

The  secant  modulus  verus  percentage  strain  at  different  confining  pressures  is 
shown  in  Figure  143.   Figure  144  shows  the  secant  modulus  at  0.1%  strain  versus 
the  confining  pressure.  At  2  psi  (0.14  kg/cm  )  pressure,  the  secant  modulus  is 
600  psi  (42  kg/cnr). 

CONSOLIDATION  TESTS 

The  ASTM  standard  method  (D2435-70)  was  followed.  The  data  from  the 
laboratory  consolidation  test  make  it  possible  to  plot  a  stress-volume  strain 
curve,  which  often  gives  useful  information  about  the  pressure  history  of  the 
soil . 
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Figure    136. 


Applied  stress  versus  average  *  strain  for 
unconfined  comcression  test. 
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Figure    137. 


Secant  modulus  versus  3  strain  for  unconfined 
compression  test. 
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Figure   140.      Modulus  of  elasticity  versus  confining 
stress  from  tri axial  compression  test. 
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Figure    143.       Secant  modulus  versus  %  strain  at  different 
confining  pressures  from  tri axial   compression  test. 


208 


6000 


A    FROm'uc"  TEST 

o   from'cu*  TFST 


1  psi  »   0.07    Jcg/ca' 


40  60  SO 

CONFINING    PRESSURE    (psi.J 


Figure  144. .     Secant  modulus  at  0.1  %  strain  versus 
confining  pressure  from  triaxial  compression  test. 
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Plots  of  coefficient  of  consolidation  (C  )  primary  compression  ratio 
(t),  and  void  ratio  (e)  versus  log  of  pressure  are  shown  in  Figure  145.   Figure 
146  shows  a  plot  of  applied  pressure  versus  %   strain.  A  plot  of  secant  modulus 
at  different  percent  strains  is  shown  in  Figure  147.   The  overconsolidation 
ratio  (OCR)  of  this  silty  clay  was  established  as  30~  hence  this  soil  was 
classified  as  an  overconsoli dated  soil. 

SOIL  CLASSIFICATION  ANO  SU8GRA0E  SEHAVIOR 

The  Atterberg  limits  and  indices  were  found,  and  grain  size  analysis 
and  hydrometer  tests  were  performed  (Lambe,  1951).  Figures  148  and  149  show 
the  grain  size  analysis  and  the  flow  curve  respectively,  the  latter  providing  the 
Plasticity  Index  of  16  for  the  soil  used.  Figure  150  shows  the  dry  density 
versus  water  content  relationship  -  the  optimum  water  content  being  3%.  The  soil 
was  mixed  with  an  over-the-optimum  water  content  (16%). 

It  is  reasonable  to  assume  that  the  plate  load  test  will  provide  a 
better  in  situ  value  of  the  modulus  of  elasticity.  The  plate  load  test  was 
conducted  in  the  sampling  area  which  did  not  undergo  as  many  preloadings  as  the 
main  experimental  area.  The  preloadings  will  increase  the^modulus.  The  adopted 
modulus  for  an  ana 


area.  The  preloadings  will  increase  the-mc 
lytical  comparison  was  850  psi  (60  kg/ air). 


From  the  consolidation  test,  the  overconsolidation  ratio  (OCR)  was 
found  to  be  30. 

From  Figure  151  ,  the  following  was  obtained:  8-10 

K0  •  2.S  (K0  •  2l) 

The  obtained  value  of  K     is  in  agreement  with  that  in  the  published  literature 
for  overconsoli dated  clays  and  compacted  clays  (Wu,  1975).     An  unacceptable  value 
of  Poisson's  ratio  is  obtained  from  the  above  K     value,  based  on  available 
relationship  for  isotropic,  homogeneous  and  elastic  material: 

(K0  -        v       ) 

0       1  -  v  B-ll 

Consequently  in  all  analytical  work  a  value  of  v  =  0,4  was  adopted.  This  value 
happens  to  agree  with  that  obtained  by  using  the  formula  KQ  =  1  -sin$,  for 
normally  consolidated  material  ($  =18  ). 
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Figure  145.     Coefficient  of  consolidation  (top),  primary  compression 
ratio  (middle),  and  void  ratio     (bottom)  at  different 
pressures  using  logarithmic  fitting. 


211 


15 


12- 


— »    9 


LU 


UJ 


/ 


3  12  16 

5  STRAIN  (in. /in.) 

1   tsf.  =    0.97   kg/cni2 
1  in  =  2.54  cm 


Figure  146.       Applied  pressure  versus  %  strain. 
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APPENDIX  C:  FINITE  ELEMENT  FORMULATION 

This  section  presents  a  brief  outline  of  the  solution  technique 
using  the  finite  element  method.  Emphasis  is  placed  on  thermal  stress  analysis 
and  transient  heat  flow  analysis.  Stiffness  matrix  derivations  are  given 
and  a  discussion  of  the  methods  used  to  solve  the  system  of  linear  equations 
and  the  system  of  partial  differential  equations  is  also  presented. 

CLASSICAL  METHOOS 

Most  classical  methods  of  plate  analysis  are  derived  from  the  equa- 
tions of  the  theory  of  elasticity.  Only  a  small  portion  of  small  deflection 
plate  problems  possess  rigorous  close-form  solutions;  as  a  result,  numerical 
methods  are  usually  applied. 

The  structures  to  be  analyzed  in  this  study  are  relatively  complex. 
The  anchored  pavement  consists  of  a  slab  with  a  non-uniform  distribution 
of  flexural  rigidity.  The  problem  can  be  modeled  mathematically  with  the 
classic  plate  on  elastic  foundation  formulation.  From  the  plate  bending 
equation 

32tji     3^m    32m  -  , 

33F1  +  2  3x#  +  3y^  "   Pz(x>^  *  P*(x'*) 


where  p  is  body  forces,  p*  is  Winkler  reactive  forces,  and  mv,  mw,  and  m 
are  internal  moments  such  that 


m  = 

X 

v3x* 

Sy27 

■y" 

-0(& 

v3xx 

.   32w> 
+  V3F} 

mxy  = 

•  -d-v) 

p32w 
3x3y 

x'  y'     xy 

C-2 

C-3 

C-4 

where  D  is  the  slab  flexural  rigidity,  v  is  Poisson's  ratio,  by  substituting 
the  internal  moment  expressions  and  keeping  D  variable,  then 

lx*  ^IT2*  +  V}   2(1"v)3^ay  °3xay  +  W  D{W 

This  governing  equation  is  not  conveniently  solved  and  may  not  possess  a 
closed-form  solution,  except  for  simple  boundary  conditions  and  slab  geometry. 
Because  the  small  deflection  theory  of  plates  ignores  second  and  higher  order 
terms,  numerical  solutions  may  be  prone  to  error.  In  addition,  this  classical 
plate  equation  with  Winkler  foundations  considers  only  vertical  reactive 
pressures  that  may  not  truly  characterize  the  behavior  of  the  system.  The 
anchored  pavement,  in  particular,  utilizes  vertical  members  which  are  restrained 
from  rotation  by  passive  earth  pressure.  This  imposes  a  rotational  constraint 
at  the  anchor  locations.  The  anchors  are  expected  to  transfer  a  portion 
of  the  applied  loads  deeper  into  the  subgrade  where  confining  pressures  are 
higher  and  as  a  result,  modulus  is  greater.  The  problem  becomes  much  more 
difficult  with  varying  subgrade  reaction  and  rotational  constraints. 
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The  finite  element  method  can  account  for  these  conditions  and 
also  will  readily  determine  stress  distributions  within  the  slab  and  the 
subgrade. 

STATIC  ANALYSIS 

The  finite  element  method  involves  the  division  of  a  continuum 
into  an  assemblage  of  discrete  continua  which  are  interconnected  at  nodal 
points  such  that  the  original  character  of  the  structure  and  its  boundary 
conditions  are  maintained.  The  formulation  of  element  stiffnesses  that  link 
applied  nodal  forces  with  nodal  displacements  through  a  relation  of  the  form 

{Fe}  *    Ke  (ue}  .  C-6 

where  (F  }  is  a  vector  of  applied  element  forces,  K   is  the  element  stiffness 
matrix,  and  {u  }  is  a  vector  of  nodal  displacements,  begins  with  the  choice 
of  suitable  functions  that  represent  nodal  displacements  in  terms  of  nodal 
coordinates.  The  displacement  functions  for  the  various  elements  used  in 
this  study  are  polynomials  and  listed  in  Table  23..  The  elements  used  are 
planar  quadrilateral  and  hexahedral  solid  elements  with  translational  and 
temperature  degrees  of  freedom. 

ANSYS  uses  a  frontal  solution  for  the  system  of  simultaneous  equa- 
tions. The  solution  proceeds  by  advancing  along  increasing  element  numbers 
rather  than  nodes  as  in  a  conventional  Gaussian  scheme  (Irons,  1970).  An 
element  with  its  associated  nodal  degrees  of  freedom  is  placed  in  core  and 
the  nodes  comprising  this  element  are  scanned.  If  the  degrees  of  freedom 
inherent  to  these  nodes  are  no  longer  associated  with  any  elements  yet  to 
be  processed,  then  those  degrees  of  freedom  are  condensed  from  the  active 
equations  in  core  by  Gaussian  elimination,  and  the  wave  front  moves  to  subse- 
quent elements.  The  solution  proceeds  until  all  elements  have  been  considered. 
This  technique  can  often  be  more  efficient  than  a  banded  equation  solver. 
ANSYS  also  has  wave  reordering  capability  that  aids  in  the  most  efficient 
formulation  of  the  structural  stiffness  matrix  (least  core  space  requirements) 
(Kohnke,  1977;  Meyer,  1973). 

The  global  formulation  of  the  problem  after  the  boundary  conditions 
are  applied  on  {u}  and  {F}  is  the  partitioned  form 

*   Kr    n_    ,  F_  C-7 

r  rr        r     r 

where  the  subscript  r  refers  to  the  applied  boundary  conditions.  -The  solution 
of  the  upper  half  of  this  equation  results  in  an  expression  for  the  unknown 
displacements 

{u}  *   -  K  ~x  Kp  {ur}  +  K  "l  F  C-8 

while  solution  of  the  lower  half  yields  reactive  forces 

{%}'  *  Kr  T{u}  ♦  Krr  {ur}  C-9 
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Strains  are  then  calculated  from  derived  strain-displacement  relations,  and 
stresses  are  determined  from  stress -strain  relations. 

TRANSIENT  HEAT  TRANSFER  ANALYSIS 

The  solution  to  the  Laplacian  field  equation 

c272u  a  u  C-10 

where  7  is  the  Laplacian  operator  and  the  dot  indicates  a  time  derivative, 

by  the  finite  element  method,  was  first  proposed  by  Zienkiewicz  and  Cheung 

(1965).  The  heat  conduction  form  in  three  dimensions  is  the  following  quasi- 
linear  parabolic  equation: 

pc3t   3x^3xj   3yl<3yJ   3z^3zJ 

where  p  is  density,  c  is  specific  heat,  k  is  thermal  conductivity,  and  T 

is  temperature.  The  finite  element  method  was  then  applied  to  heat  conductivity 

(Wilson  and  Nickel!,  1966). 

The  element  level  formulation  proceeds  with  the  statement  of  the 
principle  of  virtual  work  which  is  an  energy  balance  relation  equating  a 
virtual  change  in  internal  strain  energy  with  a  virtual  change  in  applied 
external  work.  Virtual  internal  work  is 

<SU  =  {<5G}T{QjdV  C-12 

V      v 

where 

(G}T -{&§§}  C-13 

a  temperature  gradient  vector,  and 

,QV}T  s  (Qx  Qy  Q2J  CM 

is  a  vector  of  heat  flow,  Q.  is  the  heat  flow  in  the  1-2-  direction  per  unit 
area  and  T  is  temperature  at  a  point  in  the  element.  Heat  flows  are  related 
to  temperature  gradients  through  the  conductivity  matrix,  0  such  that 

{Qv}  -  0  {G}  C-IS 

Temperature  gradients  are  related  to  nodal  temperatures  by 

{G}  -  S  {T  }  C16 


where 


{T}  =  N  T(Ta}  C-17 

8  T  .  {M}T  (M,T  (|NjT  c.w 
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and  {T}  is  a  vector  of  temperatures  within  the  element,  {T  }  is  a  nodal  tempera- 
ture vector,  and  N  is  a  matrix  of  shape  functions.  Thus 

<5tl  *   {STjT  B  T  D  B  {T  }dV  C-19 

V    e  e 

The  internal  work  associated  with  convection  surfaces  is 


6U  *  5ATQndA  C-20 

A    n 


where 


AT  =  T|s  -  T3  C-21 

with  T|  being  the  temperature  function  evaluated  at  a  point  on  the  convection 
surface  while  Tg  is  the  bulk  temperature  of  the  ambient  material.  For  a 
constant  T„, 

<5AT  =  <5T|S  C-22 

Defining  heat  flow  over  a  unit  area  as 

Qn  =  hfAT  C-23 

where  hf  is  the  film  coefficient  of  surface  heat  transfer,  yields  the  resulting 
equation 

SU  =  {STe}Thf{N|s}T{Te}dA  .   {STe}ThfT8{N|s}dA  C-24 

where  {N|  }  are  the  shape  functions  evaluated  on  the  convection  surface. 

Internal  heat  generation  rate  has  the  associated  virtual  work 

511  =    5T  q  dV  C-25 

V 


or 


where 


<SU  s     {<5Ta}T  q   {NfdV  C-26 

V       e 

The  virtual  work  due  to  a  change  in  stored  energy  is 

<5U  =     5TfdV  C-27 

V 


which  is  the  total  change  in  heat  per  unit  volume  per  unit  time.  Specific 
noted  by  c  and  density  by  p.  Thus 

5U  -  {5TJT  cn{N}{N}T{tfl}dV  C-29 
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where  {t  }  is  time  varying  temperatures. 

Consideration  of  nodal  heat  flows  gives 

<5W  ■  {6Te}T{Qe}  C-30 

where  {Q  }  is  the  nodal  heat  flow  vector  associated  with  the  element. 

Combining  the  internal  work  and  equating  this  to  external  work 
results  in  the  following: 

(  Ke  +  <|  HV  *  Ce  {ta}  »  {q<}  +  {q|}  +{qe}  c-31 

where  K   is  the  conductivity  matrix,  K^  is  conduction  matrix  associated 
with  convection,  Ca  is  the  heat  capacity  or  thermal  damping  matrix,  {Q^} 
is  internal  heat  generation,  [Qz\   is  internal  convection  heat,  and  {Q  J  Ts 
applied  nodal  heat  flow. 

The  basic  thermal  diffusion  equation 

C  .t}  +  K  {T}  =  {Q}  C-32 

in  global  form  is  solved  by  a  quadratic  temporal  integration  scheme  by  the 
following  substitutions: 

m  3  3Tt+At  '  4Tt  *  Tt-At  C-33 

Vi  2At 

{T}  =   (W  C-34 

where  t+At,  t,  and  t-At  are  current,  past  and  second  previous  past  times. 
The  solution  for  temperature  distributions  at  the  current  time  is 

W  •  ( K  +  m c >"'  <W  +  c  <4Tt  "Jt'"')  c-35 

If  K  and  C  are  temperature  dependent,  then  each  are  evaluated  at  the 
current  time  based  upon  extrapolation  from  the  previous  time. 

A  current  time  solution  requires  information  from  two  previous 
times,  thus  the  solution  procedure  is  not  self  starting.  If  the  analysis 
starts  at  time  t  s  0,  then  a  steady  state  analysis  is  performed.  If,  however, 
the  first  load  step  ends  at  time  t  >  0,  then  the  initial  two  previous  tempera- 
ture distributions  are  set  to  a  uniform  temperature  selected  by  the  user  and 
the  solution  proceeds  from  there. 
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FEDERALLY  COORDINATED  PROGRAM  (FCP)  OF  HIGHWAY 
RESEARCH  AND  DEVELOPMENT 


The  Offices  of  Research  and  Development  (R&D)  of 
the  Federal  Highway  Administration  (FHWA)  are 
responsible  for  a  broad  program  of  staff  and  contract 
research  and  development  and  a  Federal-aid 
program,  conducted  by  or  through  the  State  highway 
transportation  agencies,  that  includes  the  Highway 
Planning  and  Research  (HP&R)  program  and  the 
National  Cooperative  Highway  Research  Program 
(NCHRP)  managed  by  the  Transportation  Research 
Board.  The  FCP  is  a  carefully  selected  group  of  proj- 
ects that  uses  research  and  development  resources  to 
obtain  timely  solutions  to  urgent  national  highway 
engineering  problems.* 

The  diagonal  double  stripe  on  the  cover  of  this  report 
represents  a  highway  and  is  color-coded  to  identify 
the  FCP  category  that  the  report  falls  under.  A  red 
stripe  is  used  for  category  1,  dark  blue  for  category  2, 
light  blue  for  category  3,  brown  for  category  4,  gray 
for  category  5,  green  for  categories  6  and  7,  and  an 
orange  stripe  identifies  category  0. 

FCP  Category  Descriptions 

1.  Improved  Highway  Design  and  Operation 
for  Safety 

Safety  R&D  addresses  problems  associated  with 
the  responsibilities  of  the  FHWA  under  the 
Highway  Safety  Act  and  includes  investigation  of 
appropriate  design  standards,  roadside  hardware, 
signing,  and  physical  and  scientific  data  for  the 
formulation  of  improved  safety  regulations. 

2.  Reduction  of  Traffic  Congestion,  and 
Improved  Operational  Efficiency 

Traffic  R&D  is  concerned  with  increasing  the 
operational  efficiency  of  existing  highways  by 
advancing  technology,  by  improving  designs  for 
existing  as  well  as  new  facilities,  and  by  balancing 
the  demand-capacity  relationship  through  traffic 
management  techniques  such  as  bus  and  carpool 
preferential  treatment,  motorist  information,  and 
rerouting  of  traffic. 

3.  Environmental  Considerations  in  Highway 
Design,  Location,  Construction,  and  Opera- 
tion 

Environmental  R&D  is  directed  toward  identify- 
ing and  evaluating  highway  elements  that  affect 

*  The  complete  seven-volume  official  statement  of  the  FCP  is  available  from 
the  National  Technical  Information  Service,  Springfield,  Va.  22161.  Single 
copies  of  the  introductory  volume  are  available  without  charge  from  Program 
Analysis  (HRD-3),  Offices  of  Research  and  Development,  Federal  Highway 
Administration,  Washington,  D.C.  20590. 


the  quality  of  the  human  environment.  The  goals 
are  reduction  of  adverse  highway  and  traffic 
impacts,  and  protection  and  enhancement  of  the 
environment. 

4.  Improved  Materials  Utilization  and 
Durability 

Materials  R&D  is  concerned  with  expanding  the 
knowledge  and  technology  of  materials  properties, 
using  available  natural  materials,  improving  struc- 
tural foundation  materials,  recycling  highway 
materials,  converting  industrial  wastes  into  useful 
highway  products,  developing  extender  or 
substitute  materials  for  those  in  short  supply,  and 
developing  more  rapid  and  reliable  testing 
procedures.  The  goals  are  lower  highway  con- 
struction costs  and  extended  maintenance-free 
operation. 

5.  Improved  Design  to  Reduce  Costs,  Extend 
Life  Expectancy,  and  Insure  Structural 
Safety 

Structural  R&D  is  concerned  with  furthering  the 
latest  technological  advances  in  structural  and 
hydraulic  designs,  fabrication  processes,  and 
construction  techniques  to  provide  safe,  efficient 
highways  at  reasonable  costs. 

6.  Improved  Technology  for  Highway 
Construction 

This  category  is  concerned  with  the  research, 
development,  and  implementation  of  highway 
construction  technology  to  increase  productivity, 
reduce  energy  consumption,  conserve  dwindling 
resources,  and  reduce  costs  while  improving  the 
quality  and  methods  of  construction. 

7.  Improved  Technology  for  Highway 
Maintenance 

This  category  addresses  problems  in  preserving 
the  Nation's  highways  and  includes  activities  in 
physical  maintenance,  traffic  services,  manage- 
ment, and  equipment.  The  goal  is  to  maximize 
operational  efficiency  and  safety  to  the  traveling 
public  while  conserving  resources. 

0.  Other  New  Studies 

This  category,  not  included  in  the  seven-volume 
official  statement  of  the  FCP,  is  concerned  with 
HP&R  and  NCHRP  studies  not  specifically  related 
to  FCP  projects.  These  studies  involve  R&D 
support  of  other  FHWA  program  office  research. 
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